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GENERAL POTENTIAL THEORY OF ARBITRARY WING SECTIONS 

By T. Theodoksen and I. E. Garrick 


SUMMARY 

Thif^ report gives an exact treatment of the problem oj 
determining the 2 -dimensional potential flow around 
wing sections oj any shape. The treatment is based 
directly on the solution oj this problem as advanced by 
Theodorsen in iV. A. C. A. Technical Report No. Jftl. 
The problem condenses into the compact form of an inte- 
gral equation capable oj yieldmg numerical solutions by 
a direct process. 

An attempt has been made to analyze and coordinate 
the results of earlier studies relating to properties oj iving 
sections. The existing approximate theory oj thin wing 
sections and the Joukowsky theory with its numerous 
generalizations are reduced to special cases oj the general 
theory oj arbitrary sections^ permitting a clearer perspec- 
tive oj the entire field. The method not only permits the 
determination oj the velocity at any point oj an arbitrary 
section and the associated lift and moments, butjurnishes 
also a scheme for developing new shapes oj preassigned 
aerodynamical properties. The theory applies also to 
bodies that are not airjoils, and is oj importance in other 
branches oj physics involving potential theory. 

INTRODUCTION 

The solution of the problem of determining the 
2-dimensional potential flow of a nonviscous incom- 
pressible fluid around bodies of arbitrary shape can be 
made to depend on a theorem in conformal represen- 
tation stated by Riemann almost a century ago, 
known as the fundamental theorem of conformal rep* 
resen tation. This theorem is equivalent to the state- 
ment that it is possible to transform the region 
bounded by a simple curve into the region bounded by 
a circle in such a way that all equipotential lines and 
stream lines of the first region transform respectively 
into those of the circle. The theorem will be stated 
more precisely in the body of this report and its sig- 
nificance for wing section theory shown—suffice it at 
present to state that if the analytic transformation by 
which the one region is transformed conformally into 
the region bounded by the circle is known, the poten- 
tial field of this region is readily obtained in terms of 
the potential field of the circle. 

A number of transformations have been found by 
means of which it is possible to transform a circle into 


a contour resembling an airfoil shape. It is obviously 
true that such theoretical airfoils possess no particular 
qualities which make them superior to the types of more 
empirical origin. It was probably primarily because 
of the difficulty encountered in the inverse problem, 
viz, the problem of transforming an airfoil into a 
circle (which we shall denote as the direct process) 
that such artificial types came into existence. The 
2-dimensionftl theoretical velocity distribution, or what 
is called tlie flow pattern, is known only for some 
special symmetrical bodies and for the particular class 
of Joukowsky airfoils and their extensions, the out- 
standing investigators * being Kutta, Joukowsky, and 
von Mises. Although useful in the development of 
airfoil theory these theoretical airfoils are based solely 
on special transformations employing onh^ a small 
part of the freedom permitted in the general case. 
However, they still form the subject of numerous 
isolated investigations. 

The direct process has been used in the theory of 
thin airfoils with some success. An approximate 
theory of thin wing sections applicable only to the 
mean camber line has been developed ^ by Munk and 
Birnbaum, and extended by others. However, at- 
tempts ^ which have been made to solve the general 
case of an arbitrary airfoil by direct processes have 
resulted in intricate and practically unmanageable 
solutions. Lamb, in his “Hydrodynamics” (reference 
1, p. 77), referring to this problem as dependent upon 
the determination of the complex coefficients of a 
conformal transformation, states: “The difficulty, 
how^ever, of determining these coefficients so as to 
satisfy given boundary conditions is now^ so great as 
to render this method of very limited application. 
Indeed, the determination of the irrotational motion of 
a liquid subject to given boundary conditions is a 
problem whose exact solution can be effected by direct 
processes in only a very few cases. Most of the cases 
for wLich wx know the solution have been obtained by 
an inverse process; viz, instead of trying to find a 
value of or ^ which satisfies [the Laplacian] = 0 
or 8-nd given boundary conditions, w'e take 

some knowm solution of these differential equations 

• See bibliography given in reference 9, pp, 24, 84, aii<l .')8L 

* CL footnote 1 

’ See Appendix II of this pai)er. 
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and inquire what boundary ronditions it can be nmdo 
to satisfy.” 

In a report (reference 2) recently published by the 
National Advisory ('oiiiinittee for Aeronautics a 
eral solution employing a direct method was briefly I 
p;iven. It was shown that the ]U'oblein could be stated 
in a condensed form as an integral equation and also 
that it was possible to elTect the practical solution of 
this equation for the case of any given airfoil. A 
formula giving the velocity at any ])oint of the surface 
of an arbitrary airfoil was develo])ed. The first part 
of the present pa]>er includes the essential develop- 
ments of reference 2 and is dt'voted to a more com- 
plete and precise treatment of the method, in ])articu- 
lar with respect to the evaluation of the integral 
equation. 

In a later part of this paper, a geometric treatment 
of arbitrary airfoils, coordinating the results of earlier 
investigations, is given. S])ecial airfoil ty]>es have 
also been studied on the basis of the general method i 
and their relations to arbitrary airfoils have been 
analyzed. The solution of the inverse problem of 
creating airfoils of special types, in particular, types of ■ 
specified aerodynamical ])ro])erties, is indicated. 

It is hoped that this ])aper will serve as a step 
toward the unification and ultimate sim])lification of 
the theory of the airfoil. 

TRANSFORMATION OF AN ARBITRARY AIRFOIL INTO 
A CIRCLE 

Statement of the problem. — The problem which this 
report proposes to treat may be formulated as follows. 
Given an arbitrary airfoil^ inclined at a specified angle ; 
in a nonviscous incompressible fluid and translated 
with uniform velocity V. To determine the theoreti- 
cal 2-dimensional velocity and pressure distribution at 
all points of the surface for all orientations, and to 
investigate the properties of the field of flow surround- 
ing the airfoil. Also, to determine the important 
aerodynamical parameters of the airfoil. Of further 
interest, too, is the problem of finding shapes with 
given aerodynamical properties. 

Principles of the theory of fluid fiow.~We shall 
first briefly recall the known basic principles of the 
theory of the irrotational flow^ of a frictioriless incom- j 
l^ressihle fluid in two dimensions. A flow is termed ! 
‘‘2-dimensional” when the motion is the same in all 
planes parallel to a definite one, say xy. In this case 
the linear velocity components u and v of a fluid 
clement are functions of x, y, and i only, 

'^Jdie differential equation of the linos of flow in this 
<’use is 

V dx-u dy 0 

* Hy ml airfoil .shape, or w ing seel ion, isrotj^jlily meiui( an olorii:a‘e<l smooth shape 
rounded at the leading edire ami enilirij: in a sharp edi^e at ttie rear. All praeUcal 
airfoils are eharaeterin'd hy alai k of abrupt eha!ij:e tdi tirvaf lire exoei't for h rumided | 
nose and a small ratlins of turvainre at the (ail. 


and the (‘(piatioii of continuity is 

du _ d{- v) 

dx ' dy dx~ dy 

which sliows tlmt the above first equation is an e.xact 
differetitial. 

If Q-c is the integral, then 


M 

dy 


and = 




This function Q is called tlic stream function, and 
the lines of flow', or streamlines, are given by the ecj na- 
tion Q c, where c is in general an arbitrary function 
of time. 

Furthermore, we note that the existence of the 
stn‘nm function does not depend on whether the motion 
is irrotational or rotational. When rotational its 
vorticity is 

dv du _ d~Q ^ d^Q 
^ dx dy dx^ dy^ 


which is twice the mean angular velocity’ or “rotation” 
of the fluid element. Hence, in irrotational flow' the 
stream function has to satisfy 


dx^ dy^ 




Tlien there o.xists a velocity potential P and we have 


dP_ ^dQ^ 

dx ^ dy 

dP_ ^ dQ 


The equation of continuity is now 
djc- ^ dy- 


Etjuations (1) show' that 

dP ^^^dP d(J 
dx dx dy dy 


(U 


( 2 ) 


so that the family of curves 


P = " constant, Q = constant 


cut orthogonally at all their points of intersection. 

For steady flow's, that is, flows that do not vary 
with time, the paths of the particles coincide wuth the 
streamlines so that no fluid ])asses normal to them. 
The Bernoulli formula then holds and the total pres- 
sure head II along a streamline is a constant, that is 


ji p ;/ = // 


where 7/ is the static pressure, v the velocity, an<l p 
the density. If we denote the undisturbed velocity 
at infinity V, the quantities p' ~ ^ hy 7;, and 

% p by g, the Bernoulli formula may be ex])ressed as 



( 3 ) 
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The solutions of equations (2) and (2^), infinite in 
niiinbor, represent all possif)le tv])os of irrotational 
motion of a nonviseous ineom])ressihl(‘ fluid in two 
dimensions. For a f^iven ])rol)lem there are usually 
<‘ertain sj)ecified houndary eon(liti(ms to he satisfied 
whieh may he suflieient to fix a uni(jue solution or a 
family of solutions. Tlie problem of an airfoil movino: 
uniformly at a fixed an^le of ineideiu'e in a fluid fi(dd 
is identical with that of an airfoil fi\('d in ])osit ion and 
fluid streaming uniformly past it. Oui- ])rohlem is 
then to determine the functions P and Q so that the 
vidocity at each point of the airfoil profih' has a direc- 
tion tangential to the surface (that is, tlu' airfoil con- 
tour is itself a streamline) and so tliat at infinite dis- 
tance from th(‘ airfoil the fluid has a (‘onstant velocity 
and direction. 

The introduction of the compl(‘x variable, z=x^hj, 
simplifies the ]>rohlem of determining 1* and Any 
analytic function ie(2:) of a co]n[)lex variable z, that is, 
a function of 2 ])ossessiiu 4 a, unicjiu* derivative in a 



eacli real functions of x and y. Suppose now in the 
x]f complex plane there is traced a simple curve /(s). 
(Fi^. 1.) Each value of s along the curve defines a 
point re in the iv ])lane and/(2) maj)S into a curve/(uO 
or P{z). Because of the special properties of analytic 
I functions of a com]>lex variable, there exist certain 
■ s])C(‘ial relations between /(c) and F{z). 

The outstanding property of functions of a complex 
variable analytic in a region is the existence of a uni(iue 
derivative at every ])oint of the region. 

d?c Pm w — , 

r = 

; (\z £— z — z 

or 

du? 

This relation expresses the fact tliat any small curve 
I zz' through the point £: is transformed into a small 
curve unr' througli the point w by a magnification p 
and a rotation y; i. e,, in Figure 1 the tangent t will 
coincide in direction with T by a rotation 7 - / 3 — a. 



Fn.t lu; 2 — r^rthojrnnfil network ob(aine<! hy a eonfornml ImnsfornKition 


region of the complex plane, may be s(‘pai“ated into its 
I’eal and imaginary ]>arts as iv(z)^ w{x ^ iy) P{x, y) 
■riQix, y)y determining functions P and Q whi(di may 
represent the veh)city potential and stream function of | 
a possible fluid motion. Idnis, analytic functions of a j 
complex variable jiossess the s])C(‘ial property that tlie | 
component functions P and Q satisfy the Cauchy-'' 
Kiemann equations (eq. (1)), and each therefore also ; 
satisfies tlie e([uation of Laplace (eq. (2)). Conversely, 
any furndion P{Xfy)-\-iQ{x,y) for which and Q 
satisfy relations (1) and (2) may be written as w{xF 
iy) = w{z). The essential difliculty of the ])roblem is 
to find the ])articular function which satisfies the 
s[)(H“ial boundary-flow conditions mentioned above for 
a sjiecified airfoil. 

The metho<l of (Tinformal rejiresentation, a geomet- 
] ic application of the complex variable, is well adapted 
to this problem. The fundamental ])ropcrties of trans- 
formations of this type may be stated as follows: 
Consider a function of a complex variable z — x^riy^ 
say e (r) analytic in a given region, such that for each 
vfdue of Zy ?e(r) is uniquely defined. The function 
ir(:r) may lx* ex{)r(*ssed as irj wbei'c ^ and r\ are ^ 


This is also true for any other ]>air of corresponding 
curves through c and w, so that in general, angles 
between corresponding curves are iireservcd. In ])ar- 
ticular, a curve zz" orthogonal to zz' transforms 
into a curve vnv" orthogonal to ww\ 

It has been seen that an analytic function /(^) may 
be written ]^{x, y) -y iQ{x, y) where the curves P = con- 
stant and (2 ' form an orthogonal system. 

If then J{z) is transformed conformally into j{w) 
-- C(c v) 1- v) that is into f[w(z)]-F(z) - R(x, y) 
: iS(x, y), the curves P(x, y) = constant, Q(x, y) = con- 
stant map into the orthogonal network of curves 
R{x, y) - constant, S{Xy y) = constant. (Fig. 2.) If the 

magnification =p is zero at a point w, tlie trans- 
formation at that point is singular and ceases to be 
conformal. 

We may use the method of conformal transforma- 
tions to find the motion about a complicated boundary 
from that of a simpler boundary. Suppose w{z) is a 
function which corresponds to any definite fluid motion 
in the z ])lane, for instance, to that around a circle. 
Now if a new variable ^ is introduced and z set equal 
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to any analytic function of say then w{z) 

becomes or ir(n representing a new motion in 

the ^ plane. This new motion is, as has been seen, 
related to that in the ^ plane in siicli a way that the 
streamlines of the 2 plane are transformed by 
into tlie streamlines of the s' plane. Thus, the con- 
tour into which the circle is transformed represents 
the profile around wliich the motion H (^) exists. The 
problem of determining the flow around an airfoil is 
now reduced to finding the ])roper conformal transfor- 
mation which ma])s a curve for which the flow is known 
into the airfoil. The exi'^tfuce of such a function was 
first sliown by Kieniann. 

e shall first formulate the theorem for a simjdy 
connected region ^ bounded by a closed curve, and 
then show how it is readily a])plied to the region 
external to the closed curve. Tlie guiding thought 
leading to the tlieorern is simple. We have seen that 
an analytic function may transform a given closed 
region into another closed region. But suppose we 
are given two separate regions bounded by closed 
curves — does tliere exist an analytic transformation 
which transforms one region ('onformally into the 
other? This (piestion is answered by Kiemann's 
theorem as follows: 

Riemann’s theorem. — The interior T of any simply 
connected region (whose boundary contains more than 
one point, but we shall be concerned only with regions 
having closed boundaries, the boundary curve being 
composed of piecewise differentiable curves [Jordan 
curve], corners at which two tangents exist being per- 
mitted) can be mapped in a one-to-one conformal 
manner on the interior of the unit circle, and the 
analytic® function which consuimnates this 

transformation becomes unique when a given interior 
point 2 o of T and a direction through Zq are chosen to 
correspond, respectively, to the center of the circle and 
a given direction through it. By this transformation 
the boundary of T is transformed unupiely and con- 
tinuously into the (ire u inference of the unit circle. 

The unit circle in this theorem is, of course, only a 
convenient normalized region. For suppose the re- 
gions in the i" plan(‘ and T> in the ?c plane are 
transformed into the unit (ircie in the 2 plane by 
^=J{z) and w=-^F{z), res])cctively, then is trans- 
formed into T> by ol)tained b}^ eliminating c 

from the two transformation eejuations. 

In airfoil theory it is in the ivgion external to a closed 
curve that we are interested!. Such a region can be 
readily transformed conformally into the region in- 
ternal to a closed curve l>y an inversion. Thus, let us 
supj)ose a point 2 „ to be within a closed curve B wliose 

’ A region of the plane is simply cocuiecfed when any closed confour lyinR 

OJitirely within the region nmy he shrunk toa ptnnt w if lumt passinn out of the region. 
Cf. reference C p. where a pro<»f of the theoreiti haserl on (Jreen’.s function is 
given. 

f Attention i.s here directed to tlie fact that an analytic function is develttpalde at 
a point in a power series convergent in any circle about the point atid etdirely 
within tlte region. 


external region is F, and tlien choose a constant k 
such that for every point 2 on the boundary of T, 
Then the inversion transformation w = 
k 

will transfonn every point in the external region 

r into a point internal to a closed region lying 
entirely within B, the boundary B mapping into the 
boundary of F', the region at infinity into the region 
near Zq. We may now restate Biernann^s theorem as 
follows: 

One and ordy one analytic function i^^/(z) exists by 
means of which the region F external to a given curxT 
B in tlie plane is transformed conformally into the 
region external to a circle C in the 2 plane (center at 
2- 0 ) su(*h that the point 2 = 00 goes into the point 


^--00 and also at infinitv. 

(I2 


This function can 


be develo])ed in the external region of C in a uniformly 
convergent series with comj)lex coefheients of the form 


r - - J{ z) -2-1- ‘ ' ( 4 ) 


by means of which the radius R and also the constant 
m are completely determined. Also, the boundary B 
of F is transformed continuously and uniquely into the 
circumference of C. 

It should be noticed that tbe^ transformation ( 4 ) is 
a normalizcni form of a more general series 


” JH - Uq 'I" 4 “ 4 * + 

2 2 

and is obtained from it by a finite translation by the 
Vf'clor Uo ftinl a rotation and expansion of the entire 
field depending on the coellicient a_,. The condition 
= l is necessary and sidHcient for the fields at 
infinity to coirndde in magnitude and direction. 

The constants Ci of the transformation are functions 
of the shape of the boundary (uirve alone and our 
problem is, really, to determine the comjilex coctli- 
I dents defining a given shape. With this in view, we 
proceed tirsi to a convenient intermediate trans- 
formation. 

The transformation ^=^2'+^^ — This initial trans- 

2 

formation, although not essential to a purely mathc- 
i matical si^lution, is nevertheless very useful and 
important, as will be seen. It represents also the key 
transformation leading to Joukowsky airfoils, and is 
i the basis ol nearly all approximate theories. 

; Let us define the jmints in the f plane l)v 
using rectangular coordinates {x, y), and the points in 
the 2' plan(‘ by 2'=ac’^^‘® using polar coordinates 
(ue’^, B). Tlie constant a may conveniently be con- 
sidered unity and is added to preserve* dimensimis. 
We lun e 
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and substitnting 2 ' ^ 

we obtain f = 2a cosh (^ + iO) 


or r = 2a cosh yp cos B + 21 a sinh \p sin 6 

Since the coordinates (j, y) are given by 


X = 2a cosh yp cos 6 
y=^2a sinh yp sin 6 


(6) 


If ^^0, then = and ^ = 2a cos 0. Tliat is, if F 
and P' are corresponding points in the ^ and z' planes^ 
respectively, then as F traverses the x axis from 2a to 
— 2a, F' traverses the circle from 0 = 0 to 0 = 7 t, 
and as F retraces its path to ^ = 2a, F' completes tlie 
circle. The transformation (5) then may be seen to 
map the entire ^ plane external to the line 4a uniquely 
into the region external (or internal) to the circle of 
radius a about the origin in the z' plane. 

Let us invert equations (6) and solve for the elliptic 
coordinates yp and $. (Fig. 3.) We have 




KuiURR 3. — Transformation by eliij)tic coordinates 


cosh yp = 


X 

2a cos 6 


sinh 


rp = 


y 

2a sin 0 


and since cosh ^^ — sinh ^yp=\ 

Y-f - y Y^i 

\2a cos 0/ \2a sin 0/ 

or solving for sin^0 fwhicli can not become negative), 


2 sin^ e = p + ^ + (7) 

Similarly we obtain 

Y + { ^ Y = 1 

\2a cosh yp) \2a sinh yp ) 
or solving for sinh ^yp 

2 sinh V = - + yj p‘ + (^y (8) 

We note that the system of radial lines 0 = constant 
become con focal hyperbolas in the ^ plane. The circles 
^ = constant become ellipses in the ^ plane with major 
axis 2a cosh yp and minor axis 2a sinh \p. These orthog- 
onal systems of curves represent the potential lines and 
streamlines in the two planes. The foci of these two 
confocal systems are located at (±2a, 0). 


Equation (8) yields two values of yp for a given 
point (x, ?/), and one set of these values refers to the 
correspondence of (x, y) to the point (ae^, 0) external to 
a curve and the other set to the correspondence of 
(x, y) to the point (ac“’^, — 0) internal to another curve. 
Thus, in Figure 3, for every point external to the 
ellipse El there is a corresponding point external to the 
circle Ci, and also one internal to C) , 

The radius of curvature of the ellipse at the end of 

the major axis is p = 2a or for small valuer, of yp, 

p^2ayp'^. The leading edge is at 

2a cosh ^^2a^l + ^^^2a+ 

Now if there is given an airfoil in the plane (fig. 4), 
and it is desired to transform the airfoil profile into a 
curve as nearly circular as possible in the 2 ' plane by 
using only transformation (5), it is clear that the axes 
of coordinates should be chosen so that the airfoil 
appears as nearly ellipticjil as possible with respect to 
the chosen axes. It was seen that a focus of an 
elongated ellipse very nearly bisects the line joining 
the end of the major axis and the center of curvature 
of this phint; thus, we arrive at a convenient choice of 
origin for the airfoil as the point bisecting the line of 
length 4a, which extends from the point midway be- 
tween the leading edge and the center of curvature of 
the leading edge to a point midway between the 
center of curvature of the trailing edge and the trailing 
edge. This latter point practically coincides with the 
trailing edge. 

The curve defined by resulting in the 2 ' 

plane, and the inverse and reflected curve defined 
by are shown superposed on the f plane in 

Figure 4. The convenience and usefulness of trans- 



FitjUKK t.— Transformation of airfoil into a nearly circular contour 


formation (5) and the choice of axes of coordinates 
will become evident after our next transformation. 



• 0 ^ • 

The transformation 2 ' = ze . — Consider the trans- 
formation 2 ' = 26 ^^^^ where/( 2 ) = ^ Each exponential 

0 ^ 

term e ^ represents the uniformly convergent series 
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whore tlio oodlicionts f are (.‘oii^)Iox iium- 

hors. For f{z) convorgont at all points in a region 
external to a eertain cirele, z' has a iinicpie real abso- 
lute value in the region and its imaginary part 

is deliiiitely defined ex(‘opt foi' integral multiples of 
2tti. When = CD, = "Plip constant o, = 

Bid is then the determining fa(Tor at infinity, for the 
atleld t infinity is magnified by and rotated by the 
angle Bq. It is thus clear that if it is desired that the 
regions at infinity he identical, tliat is, z^ = z at infinity, 
the constant Cu must he z(‘ro. The constants Ci and (v 
also play im])ortant roles, as will he shown later. 

We shall now transform the closed curve ^ z' 
into the cinde z — (radius origin at center) 

hv means of the p;eneral transformation i rc’h’nuice 2) 


i27T 

-^0 


(c) 


V 


( 10 ) 


wliich leaves the fields at infinity unaltered, and we 
sliall obtain (‘X]>ressions for t!ie c'OTistants B„y and 
i/'o- The justifi<'ation of the solution will ho assured by 

the actual (‘onveiirence of since if the solution 

1 - 

exists it is unitpie. 

By definition, for the corresi)ondence of tfie houml- 
nry points, we luive 

+ v^) (10') 


Idle evaluation of tlie infinite number of constants 
as rei)resented by equations (a) and (b) can he made 
to depend upon an important single equation, whicii 
we shall obtain by eliminating these constants from 
equation (12). 

Substitution of (a) and (l>) for the coefficients of 
e(j nation (12) gives 


{6 — (f)' = " 2: 1 cos Uip' J' \l/{(p) sin Uip d^ 
0 


-sin }}(p' f cos d<^ I 

« J 

whei'e 4/ (ip) \p and (9 — (p)' represents 9 — p as a func- 
tion of <p\ and where v?' is used to distinguish tlie angle 
kept (‘onstant wliile the integrations are performed. 
The exq^ression may be readily rewritten as 


(9—<p)' = - N J' \4(<p) (sin Hip coi^ Ti(p' ~ V(p sin Uip' ) dtp 

7T 1 

1 « 

- ::: f 4^{<p) sm n(ip-<p') dip 

^ 1 0 




Also 

Consequently 


cos {2n -f 1 ) 


2 sin 


ip- p 


1 

where z - - ae.'^ 

On writing ^^^/i^(cos p \ sin p\ where 7C 
have 

1 

'A - + I - <p) np - f sin np) 

Kqunting the real and imaginary parts of this relation, I 
we obtain tlie two conjugate Fourier expansions: 

'I'-io ^ ~ II V- t sin « .p J (11) 


But 

” • / t 1 L P P 

N sin n(p - p ) =- cot - 
Then 

cos(2n+l) ^ 


(p- p) 


1 “7T 

~.r,f^(p) 




.1, . 

cos np sin 




(121 


The first integral is independent of //, while the latter 
one heroines identically zero. 

Then finally, representing p — 0 by a single quantity 
e, viz p - 9^€~^(:{p)^ wc hav^e 

2tt 


1 


From equation (11), the values of the coefficients y^”> 


Br. 


R 


and the constant \pQ are obtained as follows: 

pU 1 

„ = ~ f 4^ cos -np dp 




0 


B^ \ -J . 

= - J ^ Sin np dp 

II 7T 2, 


(«) 

(b) 


iiip'} - - ,r f 'I'M cot Mr- 

By solving for the eoefhcients in equation (12) and 
substituting these in equation (11) it may be seen that 
a similar relation to equation (13) holds for the func- 
tion 4'(p)’ 


' I'lile.ss otluTwisti slaJi'ii. 4 ' inni 0 vs ill imw he us<*i1 in Uiis rt'stririt'il .<ciisi', i. (>., ji 
lipfiTiinj: tlip t oriniliiry rtirvo ilsolf, iiti'i nnt nil points in the r' [ilanp. 


i{p')=^ ,V- /c(si’) cot ”y^d<^T — y'tA(v?)ds^ (14) 

"’T() 2 

The last term is merely the constant which is, as 
has been shown, determined by the condition of mag- 
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nification of the z and 2 ' fields at infinity. The 


[ 

corresponding integral-- S €{<p) dv? docs not appear in 
2iro 


equation (13), being zero as a necessary consequence 
of the coincidence of directions at infinity and, in 
general, if the region at infinity is rotated, is a constant 
difTerent from zero. 

Investigation of equation (13). — This equation is 
of fundamental importance. A discussion of some of 
its properties is therefore of interest. It should be 
first noted that when the function 4^{<p) is considered 
known, the equation reduces to a definite integral. 
The function ® e(v?) obtained by this evaluation is the 
^‘conjugate” function to so called because of the 

relations existing between the coefficients of the 
Fourier expansions as given by equations (11) and (12). 
For the existence of the integral it is only necessary 
that \p{(p) be piecewise continuous and difTerentiable, 
and may even have infinities which must be below 
first order. We shall, however, be interested only in 
continuous single-valued functions having a period 27 t, 
of a type which result from continuous closed curves 
with a proper choice of origin. 


If equation (13) is regarded as a definite integrak it is seen 
to be related to the well-known Poisson integral which solves 
the following boundary-value ])roblein of the circle. (Reference 
3.) Given, say for the z plane a single-valued function 
for points on the circiinifcrcnce of a circle w = Re^'^ (center at 
origin), then the single-valued continuous potential function 
M(r,ff) in the external region 2 = re** of the circle which assumes 
the values u {R, r) on the circumference is given liy 


1 

uir.a) S u(R,t) 




\-r^~2Rr cos (o- — t) 


and similarly for the conjugate function 




-Ri 


0 


R-^r^ — 2Rr cos (<t--t) 


Tlicsc may be written as a single equation 


' integral equation whose process of solution becomes 
more intricate. It would be surprising, indeed, if 
anything less than a functional or integral equation 
were involved in the solution of the general problem 
stated. The evaluation of the solution of equation (13) 
is readily accomplished by a powerful method of suc- 
cessive approximations. It will be seen that the 
nearness of the curve to a circle is very signifi- 

cant, and in practice, for airfoil shapes, one or at most 
two steps in the process is found to be sufficient for 
great accuracy. 

Tlic quantities ^ and e considered as functions of 
have been denoted by ^(<^) and €(^), respectively. 
When these quantities are thought of as functions of 6 
they shall bo written as ^(0) and i{6)f respectively. 

Then, by definition 

^(d)=\p[<f{d)] (lo) 

Since — we have 

&M ■=(p-e{ip) 1 
<p{$) =$-\- i[6) i 

We arc seeking then two functions, and 

conjugate in the sense that their Fourier series expan- 
sions are given by (11) and (12), such that yl/[<p{6)] " 
where ^{6) is a known single-valued function of 
period 

Integrating equation (13) by parts, we have 


. 27T 

f log sin 
0 


d\^(<p) 
2 d<p 


(13') 


Tlie term log sin is real only in the range ip = <p' to 

v? = 27T-rv?'> but we ma}^ use the interval 0 to 27 t for ip 
with the understanding that only the real part of the 
logarithm is retained. 

Let us write down the following identity: 


u(r,a) + U'(r,ff) 

wliere the value /(z) at a point of the external region z=re*^ is 
expressed in terms of the known values /(it') along the circum- 
ference w=Re'^, In i)articular, we may note that at the 

boundary itself, since 1 ^ 77 ^^ = cot ^ \ have 

u(R,a) +iv(.R,<,)=-:rf I cot dr. 

-IQ .! 

which is a special form of equations (13) and (14). 

The quantity \p is immediately given as a function 
of 6 when a particular closed curve is preassigned, and 
this is our starting point in the direct process of trans- 
forming from airfoil to circle. We desire, then, to find 
the quantity \p as & function of p from equation (13), 
and this equation is no longer a definite integral but an 

‘ Thia function will be talle-i "cotjfDnnal angular diatortion” fumtion, for reasons 
CM ident later. 


log sin ^ =Iog sin ^ 


. e-e' 


. ((9-r6i)-(l9+6i)' - € 2 )" (<^+< 2 )' 

sin j- — sm 

+ log + log T-(-e-+e,)-(e+i,)' 

Sin —7; — sin " 


+ - log ■ 


(6 + t*) “ (6 + ' 

^ 2 


(17) 


. {6 F €n) — 


. {9 + €) “ 


sm ' 


-i- log 


Sill .) 


where in the last term we recall that e + i(9) — ; and 

where it may be noted that each denominator is the 
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niinierator of tlie prccedinj>: term. The symbols 
(4’ = 1, 2, n) represent functions of 0, which tlius 
far are arbitrary.^ 

Since by equation (15) ^(0) = we have for 

corresponding elements and d</? 

d'4^(e) , 

Then multiplying the left side of equation (17) by 
X right side ]>y ~ inte- 

giating over tlie ])eriod 0 to 2t we obtain 


e[^{0')]=€{e')=-- f log sin^^ 

7T 0 2 


dd 


sin 


d^(e) 




,, ■ (g -| -tW) - (g + e(9))' 

+ 3 y log T^z^-rArn^r-^ (18) 




. (e^i„)~(e + t„y 

Sin - -- 


where /: = 1, 2, . . 7i. 

We now' choose the arbitrary functions e*(0') so that 


and 


h{d')=- f log sin 

^0 


. (^+. n -.)-(0 T - dv ^( 0 ) 


d^ 


d^ (19) 


where 1, 2, . . w. 

Equation (18) may then be w'ritten 

e(^') =^ + <1 + («2” «i) . . ■ +(€„-- c„_i) + (« — €„) (20) 

or €(^') = Xi+\ 2 + , . . X„+X 

wdiere X*(0') ^ €*_i and is in fact the term of 

equation (18). The last term we denote by X. 

From equation (19) we see that the function t*(^?') is 
obtained by a knowledge of the preceding function 
€*_i(0O* For convenience in the evaluation of these 
functions, say 

m-if log .in » - j, 

w'e introduce a new variable v?* defined by 

^*(0) = 0+ «a(0) (^^1,2,. . .,71.) 

Then 




= 1 j log sin 

^ 0 ^ dipt 


From the definition of pt as 

Pt(e)^e+h(e) 


( 21 ) 


* The symbol rejireseiits ami is used to denote the same function 

of S' that is of 6. The variattles 6 and are reRarded as in«lei>endent of each 

other 


w'c may also define the symbol €jt(v?*) by 


wdiere 


^(<^jt) = ^k~€k(<Pk) 


^kW =€dy»ifc(0)] 


It is important to note that the symbols t*, €*, e** 
denote the same quantity considered, how'cver, as a 
function of 6, v?a, ^k-i, respectively. 

The quantities (e* — €/t~i) in equation (20) rapidly 
approach zero for wdde classes of initial curves ^^(0), 
i. e., ^'[0 (<^a) 1 very nearly equals \k[6(<pn^i)] for even 
small ^*’s. The process of solution of our problem is 
then one of obtaining successively the functions 
»A[^(v’i)]r ^[^(v? 2 )l, .... ^'lOifPn)] where 4^[6{(pr,)\ and 
«n[^(¥7/j)l become more and more ‘‘conjugate.'' The 
process of obtaining the successive conjugates in prac- 
tice is explained in a later paragraph. We first pause 
to state the conditions w'hich the functions <pt are sub- 
ject to, necessary for a one-to-one correspondence of 
the boundary points, and for a one-to-one corre- 
spondence of points of the external regions, i. e., the 
conditions wdiich are necessary in order that the 
transformations be conformal. 

In order that the correspondence betw'een boundary 
points of the circle in the 2 plane and boundary points 
of the contour in the z' plane be one-to-one, it is 
necessary that OM be a monotonic increasing function 
of its argument. This statement requires a word of 
explanation. We consider only values of the angles 
between 0 and 2tt. For a point of the circle boundary, 
that is, for one value of v? there can be only one value 
of d, i. e., 0{<f) is ahvays single valued. How^ever, (p{6), 
in general, does not need to be, as for example, by a 
poor choice of origin it may be many valued, a radius 
vector from the origin intersecting the boundary more 
than once; but since we have already postulated that 
^(0) is single valued this case can not occur, and ip{6) 
is also single valued. If w'e decide on a definite direc- 

d^ 

tion of rotation, then the inequality ^0 expresses 

the statement that as the radius vector from the origin 
sw'ceps over the boundary of the circle C, the radius 
vector in the z' plane sw'eeps over the boundary of B 
and never retraces its path. 

The inequality 

uv? dv? 

corresponds to 

d^M 

d^ 

Also, the condition 


dip 

dJ 




corresponds to 
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Multiplyinfr t>y \vc get 

This relation is shown in Figure 5 as a rec'tangular 



iiyperhola. We may notice then that the moriotonic 

behavior of ,f(d) and 6{ip) requires that remain on 

the lower branch of tb(‘ fn’IJerbola, i. e., 

_ de 


( 22 ) 


It will be seen later that the limiting valm^s 


de((^) ^ ^ dtfv?) _ 
dv? ’ dv? 


a> c„ 

w 


iW 

(I 2 




0 for 


' /i or since 


correspond to points of infinite velocity and of zero 
velo(“ity, resi>ectivelyj arising from sharp corners in the 
original curve. I 

The condition foi’ a one-to-one conformal corre- 
s])ondence between i)oints of the external region of the 
circle and of the external region of the ('ontoui\in the 
z' plaiu' may ]>e given (reference o, p. 98 and reference ; 
(>, Part II) as follows: There must be a one-to-one j 
boundary ])oint c(»rrespondence and the derivative of | 


the analytic function 2 ' “ given by equation (10) 
must not vanish in the region. That is, WTiting g{z) 

for ^ we have 


the integral transcendental function does not vanish 
in tile entire plane, the condition is equivalent to 

1' 'I'lu* values of tlit* upper lintiich of the liyperliohi arise when the reKiori iiitertial 
to ttie curve (if is traiisf(»rnu*(l into the exterrUil reeion of a circle, t)ut, may also 
Oiere l>e avoiilcti hy (leriniiii.’ \ B insiea'l of ip-8 
407C.N :\4 


By equation (10') we have on the boundary of the 
(‘ircle, g{Re^^) — \pQ— it, and 

, d[ iA(v?) -uM] 
d<p 


dz 


Re^^ 


dtj ip) _ -d^M 
d(p dip 


the first term on the right-hand side being real and the 
last term a ])ure imaginary. We have alrea<!y jiostu- 
lated the condition 

<- de 

^ - < 1 
i\\p 


a one-to-one boundary ])oint corre- 

dg{z) _ 


as necessary for 

spondence. Now by writing and 2 


d2 


R{^fV) 1 note that 


de(y>) 

dv? 


the houndarv 


values of a harmonic fuiK'tion and therefore this 

function assumes its maximum and minimum values 
on the boundary of the cin*le itself. (Reference 3, ]). 


22;i.) Ilonce 


externah region, 


e., 


never become — 1 in the 
(uin never vanish in this 


can 

(^ 2 ' 

(I2 

legion. 

At each step in the ])rocess of obtaining the succes- 
sive conjugates we desire to maintain a one-to-one 
corres])ondence between 6 and (p^, i. e., the functions 
6{(pt) and fpniO) should be nionotonic increasing and are 
hence subje<‘t to a restriction similar to equation (22), 
viz, 

(t 


dipk 


( 22 ') 


The process may be summed u]) as follows: We con- 
sider the function ^(0) as known, where \p{d) is the 
functional relation between \p and B defining a <dosed 
curve The conjugate of ^{6) with respect to 0 

is ei(0). We form the variable (pi^d-rt^{9) and also 
the function ^^{^{ipi)]. The conjugate of tA[^(<^i)l with 
respc(‘t to tpi is (* 2 (^ 1 ) which exjiressed as a function of 
$ is €j(0). We form the variable (p 2 = 0 ^ e,{d) and the 
function The conjugate of 4 i[d{<p 2 )] is €* 3 (v?:>), 

which as a function of 6 is 6j(0), etc. The graphical 
criterion for convergence is, of course, reached when 
the function ^[^(v?«)] is no longer altered by the 
jirocess. The following figures illustrate the method 
and exhibit vividly the rapidity of convergence. The 
numerical calculations of the various conjugates are 
obtained from formula I of the appendix. 

In Figure t), the }p(d) curve represents a circle re- 
ferred to an origin wliich bisects a radius (obtained 
from an extremely thick Joukowsky airfoil) (seep. 200) 
and has numerical voilues apj)roximately fiv^e times 
greater than occur for common airfoils. The \{^{<p) 
curv^e is known independently and is represented by 
the dashed curvT. The process of going from ^(^) to 
}p{(p) assuming as unknown is as follows: The 

function the conjugate fuiudion of ^(0), is found. 
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The quantity is then plotted against the new variable is drawn at P\ This process yields the function h(B) . 
v?i — 0 + €i(0) (i. e., each point of 4/{B) is displaced hori- The quantity ^ is now plotted against the new variable 
zontally a distance l,) and yields the curve ip[9(¥.,)]. + (i. c., each point of vt(9) is displaced hori- 

(Likewisc, i,(S) IS plotted against yielding zontally a distance h) giving the function i[e{^ 2 )\- 



Argument (Q^ ^ radians) 

Figuke G.— The iirttcess of obtaining successive conjugates 

The function is now determined as the conjugate | This curve is showm with small circles and coincides 

function of f\e{<p{)]. This function expressed as a with Further application of the process can 

function of eis t*2[i^i(0)] It is plotted as follows; I yield no change in this curve. It may be remarked 



Figure 7. — Process applied to transforming a square into a circle 

At a point F of t*2 (^,) and ^of edvi) corresponding to here that for nearly all airfoils used in practice one 

a. dehnite value of v'l one finds the value of $ which step in the process is sufficient for very accurate results 

corresponds to by a horizontal line through Q meet- As another example we shall show how a square 

ing «,(«) m Q ; for this value of 9, the quantity «2 at P (origin at center) is transformed into a circle by the 
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method. In Figure 7 the ^{6) curve is shown, and in 
Figure 8 it is reproduced for one octant. The value 
is ^(^) = log sec Q. The function is shown 

dashed; the function is shown with small 

crosses; and is shown with small circles. The 

solution is represented by the curve with small 
triangles and is obtained independently by the known 
transformation (reference 3, p. 375) which transforms 
the external region of a square into the external region 
of the unit circle, as follows: 


^ /^T 




\ + r 


where denotes a power series. Comparing this 

with equation (10), we find that \p{<p) except for the 
constant is given as the real part of log 


evaluated for 2 =e*^, and that e{<p) is given as the 
negative of the imaginary part. It ma}" be observed 
in Figure 8 that the function ^[Bi^p^)] very nearly 



equals ^ (</?). The functions f(<^) and «(0) are shown in 
Figure 7 (a); we may note that at which corre- 
sponds to a corner of the square, also, 

d€ 

— = C£> . 

dd 


’I Hecituse of the syriiinetry involved only the intcrvfil 0 to ^ need be uscfl. 'I'he 
iritepral in the appendix can he trefited as 

r 

^ '4 

= — r .Z* ^{^)|fot 2(^p—•p*)~('o^ 2(v5+v’^)]dv' 

T 0 


It may be remarked that the rapidity of convergence 
is influenced by certain factors. It is noticeably af- 
fected by the initial choice of €;,(0). The choice 
€o(0) "^0 implies that B and ip arc considered to be very 
nearly equal, i. e., that represents a nearly cir- 

cular curve. The initial transformation given b}" 
equation (5) and the choice of axes and origin were 
adapted for the purpose of obtaining a nearly circular 



KKiURe 9,— Translation by the distance OM 


curve for airfoil shapes.' If we should be concerned 
with other classes of contours, more appropriate 
initial transformations can be developed. If, how- 
ever, for a curve tlie quantity e ^fp~B has large 

values, either because of a poor initial transformation 
or because of an unfavorable choice of origin, it may 
occur that the choice eo(0)=O will yield a function 


€i(ipi) for which may exceed unity at some points, 

u^i 

thus violating condition (22^). Such slopes can be 
replaced by slopes less than unity, the resulting /wric- 
tion chosen as fo(^) and the process continued as 
before. Indeed, the closer the choice of the function 
€o(^) is to the final solution e(0), the more rapid is the 
convergence. The case of the square illustrates that 
even the relatively poor choice eo(0) =0 does not appre- 
ciably defer the convergence. 

The translation ^i = 2 + Ci. — Let us divert our 
attention momentarily to another transformation 
which will prove useful. We recall tliat the initial 
transformation (eq. (5)) applied to an airfoil in the ^ 
plane gives a curve B in the z' plane shown schemati- 
cally in Figure 9(a). Equation (10) transforms this 
curve into a circle C about the origin 0 as center and 
yields in fact small values of the quantity <p — B. We 
arc, however, in a position to introduce a convenient 
transformation, namely, to translate the circle C into 
a most favorable position with respect to the curve B 
(or vice versa). These qualitative remarks admit of a 
mathematical formulation. It is clear that if the 
cur\'e B itself happens to be a circle the vector by 
which the circle C should be translated is exactly the 
distance between centers. It is readily sho^\^l that 


The first step in Uie process is now to define = and form the function 

^p(v>o)]. The conjugate function of i‘[e{<f>vi)] is which expressed as a function 

of is etc. 

See p. 200. 
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then equation (10) should contain no constant term. 
We have 


z' = ze 
== z 


“ c 
i F" 


(i+3+ . )» 


= 2(l + ^‘ + ^^ + 


wiiere ' 


■ ) 

k2 = Ci-\- \ 


kj — C3 r Ciifi (- 


( 10 ) 


1 + -^ + 


■) 


X 


(10a) 


It is tlius apparent that if equation (10) contains no 
first harmonic term, i. e., if 

C[ = Ai + IBi = S \l/€'‘^d<p = 0 , 

^ 0 

the transformation is ol)tained in the so-called normal 
form 

^' = ^.^-r+^2+ ■ • • (-^3) 

-1 -I 

This translation can be effected either by substituting 
a new variable Zi = ^ f Ci, or a new variable = 2' — Ci. 



This latter substitution will be more convenient at 
this time. Writing 


2/ and z'^ae'^'^^^ 

we have 

The variables and di, can be expressed in terms of 
6, 7, and 5. In Figure 9(b), P is a point on the P 

“ These constants can be obtait-cJ in a res iirsinn form. See footnote 16. 


curve, i. e., OP — ae'^y PQ represents the translation 
vector Cl -- (W+^\ OQ is and angle POQ is 

denoted by Then by the law^ of cosines 

:= ^2^ 4- £^7 - 2e^e^ cos {6 - d) (aj 


and by the law of sines 


sin n = 


sin (B — 6) 


or 


&i^d^ ^ ^ ^4- tan ^ 


€2 * sin (^ — 6) 

1 — cos {B — b} 


(b) 


In Figure 10 are shown the 4'(B) and e{B) curves for the 
Clark Y airfoil (shown in fig. 4) and the i^i(0i) and 
ci(d]) curves which result when the origin is moved 
from 0 to M. It ma}^ be noted that <i(&i) is indeed 
considerably smaller than €{B). It is obtained from 


(if — 0])' = — S ypi (<^) cot ^ 


2 


and the constant i;^ given by 



2x 

f ypi{ip)dyp 

0 


The combined transformations. — It will be useful to 
combine the various transformations into one. We 
obtain from eejuations (5) and (10) an expression as 
follows: 


r ^ 2a cosh I log ^ 4- 2 ) 

y a ^ z J 


(24) 


or we can also obtain a power series development in 2 


c.+^+“'+!i+>+ 


where ‘ 


On = kr. 


2^ 


(25) 


The constants k^ mny be obtained in a convenient 
recursion form as 


ky — Cl 

2ki ^ A'lCi - 2 c2 

dA*3 - k^Ci ^ 2kxC2 4- 3c3 

4A*4 “ ■ k[\C\ 4" 2A”2C2 4’ dA*iC3 "t 4 C 4 


The constants have the same form as A'„ but with 
each c, re])laced by ~Cf (and Ao^ 1). It will be re- 


The foiislant is in\ ariant to change of origin. (.See p, 200 ) It siionbi be 
remarked that the Iranslathuj f>y tfie vector Ci is only a matter of conveiiieiKx^ and 
is esi>ediilly useful for very irregular .‘jliapes. For a study of the properties of airfoil 
shapes we shall iisi,':o;il> the oriiimil «(v?t curve. (Fig. 10(a).) 
i* By equations (,ii and (10) we have 




2 ? 
I 


-1' - 


The constant tin is thus the c(xdlicient of — in the expansion of c * an<itiie ( onstant 


hn the coetficiejit of in the expansion of e 
see Sndthsonian .Matheniatical Formula’ and 



iz'. 


For the recursion form for kn 


Tables of F.lliptic Functions, p, 120, 


268 




GENERAL POTENTIAL THEORY OF ARBITRARY WING SECTIONS 


called that the values of are ^iven by the coefficients 
of the Fourier expansion of ^(<^) as 


r 1 

7 ^= f where 11 = 

7T 0 

and 

1 2^ 

o'- 

^7T Q 

The first few terms of c(| nation (25) are then as 
follows: 


^ , , C 2 + 2 , C 3 + C 2 C, t) , 

f = 2 + Cl H 1 -—2 ^ 


(250 


By writing Zy 2 + Cl, equation (25) is cast into the 
normal form 

=" 2] + 0 + -_ 2 F • • • (“h) 

i-X -ij 


The constants 6„ may he evaluated directly in terms 
of (In or may he obtained merelv hv replacing y^M hy 
(ip) in the foregoing values for a„. 

The series given hy equations (25) and (2()) may he 
inverted and 2 or 2 i developed as a power series in 
Then 


/ % fh 

and 


CL<t ~l“ d,\C\ 


a,c ,^ ') 2«2Ci + a, -f- 0,“ 


. . . (27) 




62 _ ^3 


(28) 


The various transformations have been performed 
for the purpose of transforming the flow pattern of a 



Fir.irKK 11.- SlreHiiiliiies at nul. circle wilh zero circulation (shown hy the fill! 
lines) </= — I'sinh m sin sP = constant 

circle into the flow pattern of an airfoil. We are thus 
led immediately to the well-known problem of deter- 
mining the most general type of irrotational flow 
around a circle satisfying certain specified boundary 
conditions. 

The flow about a circle, — The boundary conditions 
to he satisfied are: The circle must be a streamline of 
flow and, at infinity, the velocity must have a given 
magnitude and direction. Let us choose the ^ axis as 
corre'sponding to the direction of the velocity at 


infinity. Then the problem stated is equivalent to 
that of an infinite circular cylinder moving parallel to 
the f axis wdth velocity V^in a fluid at rest at infinity. 

The general complex flow potential for a circle of 
radius B, and velocity" at infinity V parallel to the x 
axis is 

u,(.)=-V'(a+f)-glog| (29) 

where F is a real constant parameter, known as the 



F'iui.'rf. 12 — Streamlines ahnut circle fyr V=0 Q= — hn = constant 

circulation. It is defined as c^c,d.9 along any closed 
curve inclosing the cylinder, r, being the velocity 
along the tangent at each point. 

Writing z = Re^^^‘^ and w = P equation (29) be- 
comes 

w= "V CQ^\\(fi-\-iip) — ^(pL + i<p) (29') 

r 

or P = — V cosh M cos ^ v? 

Ztt 

Q = —V sinh pi sin <p — :^M 
Jtt 

For the velocity components, we have 



In Figures 11 and 12 are shown the streamlines for 
the cases F = 0, and V"=0, respectively. The cylinder 
experiences no resultant force in these cases since all 
streamlines are symmetrical wdth respect to it. 

The stagnation points, that is, points for w^hich u 

and V are both zero, are obtained as the roots of 4^ = 0. 

d2 

This equation has tw^o roots. 

ir ± vr67r"i?n^2_ p 2 
^ “4,ry 

and we may distinguish diflerent types of flow accord- 
ing as the discriminant F^ is positive, zero, 

or negative. We recall here that a conformal trans- 
formation w=f(z) ceases to be conformal at points 

where ^ vanishes, and at a stagnation point the flow 

divides and the streamline possesses a singularity. 

1’ Reference 4, p, or reference r>, p, 118. The log term mast be adfled because 
the region ouCside the infiniic cylinder (the point at infinity excluded) is doubly 
connected and Uierefore we must indinle the possibility of cyclic motion. 
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The different types of flow that result according 
as the parameter ^ are represented in 

Figure 13. In the first case (fig. 13 (a)), which will not 
interest us later, the stagnation point occurs as a 
double point in the fluid on the t? axis, and all fluid 
wdthin this streamline circulates in closed orbits around 
the circle, while the rest of the fluid passes downstream. 
In the second case (fig. 13 (b)), the stagnation points 
are together at <S on the circle and in the third 
case (fig. 13 (c)) they are symmetrically located on the 
circle. We have noted then that as Y increases from 
0 to At^RV the stagnation points move downward on 

the circle Re^"^ from the $ 
axis toward the rj axis. 
Upon further increase in 
r they leave the circle and 
are located on the rj axis in 
the fluid. 

Conversely, it is clear 
that the position of the 
stagnation points can de- 
termine the circulation F. 
This fact will be shown to 
be significant for wing- 
section theor}^ At pres- 
ent, we note that when 
both r and 0 a marked 
dissymmetry exists in the 
streamlines with respect to 
the circle. They are sym- 
metrical about the 7} axis 
but are not symmetrical 
about the ^ axis. Since 
they are closer together on 
the upper side of the circle 
than on the lower side, a 
resultant force exists per- 
pendicular to the motion. 

We shall now combine 
the transformation (27) 
and the flow formula for 
the circle equation (29) and obtain the general complex 
flow potential giving the 2-dimensional irrotational flow 
about an airfoil shape, and indeed, about &ny closed 
curve for which the Kiemann theorem applies. 

The flow around the airfoil. — In Figure 14 are 
given, in a convenient way, the different complex 
planes and transformations used thus far. The com- 
plex flow potential in the z plane for a circle of radius 
R origin at the center has been given as 

w(2 )= -f(s + j)-|^log2 (29) 

where F, the velocity at infinity, is in the direction of 
the negative ^ a.xis. Let us introduce a parameter to 


permit of a change in the direction of flow at infinity 
by the angle a which will be designated angle of attack 
and defined by the direction of flow at infinity with 
respect to a fixed axis on the body, in this case the 
axis ^ = 0. This flow^ is obtained simply by writing 
for z in equation (29) and represents a rotation of 



the entire flow field about the circle by angle a. We 
have 

«;(2)=-F(2e'“ + ye-'“)-|^log2 (31) 


--V.'-(l-|.--)-£ (32) 

Since a conformal transformation maps streamlines 
and potential lines into streamlines and potential lines, 



Figure 13. — Streamlines about circle 


(from Lagaliy— Handbuch der Physik 
nd. VII] Q = V sinh ^ sin = 

stant (a) n > (b) ri= 

(c) n<I6ir>ftU'i 
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we may obtain the complex flow potentials in the 
various planes by substitutions. For the flow about 
the circle in the Zi plane, 2 is replaced by 2 i — Ci 


W(2|) = - V 1^(2, - C,)«'“ + 


(i.-f.) 




c.) (31') 


(32') 


For the flow about the B curve in the z' plane, z is 
replaced by 2 ( 2 ') (the inverse of eq. (10a)) and for the 
flow about the airfoil in the ^ plane 2 is replaced by 
2 (^) from o([ nation (27) 

ir(f)=-Fl 2 (f)«'» + ^|^«--] log 2(0 (33) 


an' 



ir -] 

dzCf) 

(If 

1 iz(nv) 

' 2^2(f)]J 

d(f) 


The flow fields at infinity for all these transformations 
haA e been made to coincide in magnitude and direction. 

At this point attention is directed to two im])ortant 
facts. P'irst, in the previous analysis the original 
closed curA e may differ from an airfoil shape. The 
formulas, when convergent, are applicable to any 
closed curve satisfying the general requirements of 
the Kieinann theorem. Howe\"cr, the peculiar ease of 
numerical evaluations for streamline shapes is note- 
worthy and significant. The second important fact is 
that the parameter V which as yet is completely unde- 
termined is readily determined for airfoils and to a 
discussion of this statement the next section is devoted. 
It will be seen that airfoils may be regarded as fixing 
their own circulation. 

Kutta-Joukowsky method for fixing the circula- 
tion. — All contours used in practice as airfoil profiles 
possess the common property of terminating in either 
a cusp or sharp corner at the trailing edge (a point of 
two tangents). Upon transforming the circle into an 

I d2 

airfoil by we shall find that 

the trailing edge if the tail is perfectly sharp (or veiy 
large if the tail is almost sharp). This implies that 


is infinite at 


the numerical value of the A^elocity 


dw 

dz 


dr 


= \v\ is 

infinite (or extremely large) provided the factor |^| 

is not zero at the tail. There is but one value of the 
circulation that avoids infinite velocities or gradients 
of pressure at the tail and this fact gives a practical 
basis for fixing the circulation. 

The concept of the ideal fluid in irrotational poten- 
tial flow implies no dissipation of energy, however large 
the A^elocity at any point. The circulation being a 
measure of the energy in a fluid is unaltered and inde 
pendent of time. In particular, if the circulation is 
zero to begin with, it can never be different from zero. 


However, since all real fluids ha\^e viscosity, a better 
physical concept of the ideal fluid is to endow the 
fluid with infinitesimal viscosity so that there is then 
no dissipation of energy for finite A^elocities and pres- 
sure gradients, but for infinite velocities, energy losses 
would result. MoreoA^er, by Bernoulli’s principle the 
pressure would become infinitely negative, whereas a 
real fluid can not sustain absolute negatiA^e pressures 
and the assumption of incompressibility becomes in- 
valid long before this condition is reached. It should 
then be postulated that nowhere in the ideal fluid from 
the physical concept should the A’^elocity become 

infinite. It is clear that the factor I must then be 

zero at the trailing edge in order to avoid infinite 
velocities. It is then precisely the sharpness of the 
trailing edge which furnishes us the following basis for 
fixing the circulation. 

It AAill be recalled that the equation 0 deter- 
mines two stagnation points symmetrically located on 
the circle, the position of which varies Avith the A^alue 
of the circulation and conversely the position of a 
stagnation point determines the circulation. In this 
paper the x axis of the airfoil has been chosen so that 
the negative end passes through the trailing 

edge. From the calculation of — S (by eq. (13)) 
the A'alue of <p corresponding to any value of 0 is deter- 
mined as V? = 0 -h c, in particular at ^ = ir, ^ = x + /3, where 
/3 is the value of c at the tail and for a given airfoil is a 
geometric constant (although numerically it varies 
with the choice of axes). This angle is of consid- 
erable significance and for good reasons is called the 
angle of zero lift. The substance of the foregoing 
discussion indicates that the point 2 = 
is a stagnation point on the circle. Then for this value 
of 2 , we have by equation (32) 


or 


= - Ve^‘ 




=0 

2wz 


dw 

d^ 


-4tRVI 


2i 


= 4irRV sin (a + /3) 


(35) 


This A'alue of the circulation is then sufficient to 
make the trailing edge a stagnation point for any value 
of a. The airfoil may be considered to equip itself 
Avith that amount of circulation which enables the 
fluid to flow past the airfoil AAuth a minimum energy 
loss, just as electricity flowing in a flat plate will dis- 
tribute itself so that the heat loss is a minimum. The 
final justification for the Kutta assumption is not only 
its plausibility, but also the comparatively good agree- 
ment Avith experimental results. Figure 15 (b) shows 
the streamlines around an airfoil for a flow satisfying 
the Kutta condition, and Figures 15 (a) and 15 (c) illus- 
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trate cases for which the circulation is respectively too 
small and too large, the stagnation point being then on 
the upper and lower surfaces, re- 
spoctively. For these latter cases, 
(<^) the complete flow’ is determinable 
only if, together w'ith the angle of 
— attack, the circulation or a stag- 

(b) nation point is specified. 

Velocity at the surface. — The 
flow' formulas for the entire field 

(c) are now’ uniquely determined by 
substituting the value of V in equa- 

„ ,, , , ... tions (33) and (34). We are, how - 

FHiTRE l.'j.— (aj How With ; ' ’ 

circulation smaller than f(»r CVCr, HI a positiou tO obtaiu UlUcll 

Kutta condition; (h) flow ^hiipler and more convenient re- 

satisfying Kuttacon<lnion; » 

(cl flow with circuiaiicn latioiis for the boundary curves 
greater than tor Kutta thej„ge|^.es Indeed, wc are chicflv 

condition ' 

interested in the velocity at the 
surface of the airfoil, which velocity is tangential to 
the surface, since the airfoil contour is a streamline of 
f^ow^ The numerical value of the velocity at the 
surface of the airfoil is 


We shall evaluate each of 
equations (32 and (35) 

dw 


d{\p~u 


dtp (U .di/' .(1^ 

z'( Ae " d^ K\e \ ^ s'( *(19 \ 

ihp " ' ^ V y ‘dj/ 

de ' de 


'-cin’ 


|dw \dw 

\dz 

ds^j j 

|djr dz 

jdz' 

■ df 1 j 

lese factors 

in turn. From 

^ iAirRV 

sin (a 

±1) 

) 

2ttz 


Re‘*, and 



)-2i'Ve-'^ 

sin(a 

+ /3) 

' Ha+f) ^ -j.- 2i 

sin (a 

+ ^)1 


By c(] nation (5) 

s' = r' + and at the boundarv z' or 

2u cosh(^4- i0) 

tit - w , , .^ d{}P^ie) 

= 2u sirih(^+ i0) — " 

2 sinh(^ + 

idr ^ 

Then (sinlF^ cos^0^ cosh'V 

(sinhV T sin^0) 

and 2c“^-^/sinh^^ + sin^0) ( 

Then finally 


dw 

’ d^f; 

dz ! |dz' 

df 

’ dz 

dz' n df 


^ - 2 i W'’^[sin(a + s^^) + sin (a 4 - / 3 )] 
id 

1^1 = 2 V[sin (a + if ) — sin (a + /3)] 

\(xZ j 

In general, for arbitrary I' we find that 

|tl^- ol- • / I \ 1 1' 

k^-.l sm (« + ■?,»+ 27 ft 


To evaluate jjp: we start w ith relation (10) 


At the boundary surface 

2' = wdiere t^ip-B and 

d(4.-U)\ 

dz zV " " d^ ) 


l ^sin (a + V?) + sin (a + /5))(^ 1 -r 
== -- ^ ^ ^ - (39) 

y (sinh=^,;- + sin^9)(l +(J'^) ) 

In this formula the circulation is given by equation 
(35). In general, for an arbitrary value of F (see 
equation (3(F)), the equation retains its form and is 
given by 

F[sin(« + ^) + ^>]( ' " 

r =- . - ( 40) 

y (sinA V d sin‘^9)( 1 + ( ||^) ) 

For the spe<‘ial case r=^0, we get 
rsin(a+,5)(l 


(sinh"tA + sin'»)( 1 + (d^) ) 


Equation (40) is a general result giving the velocity 
at any point of the surface of an arbitrary airfoil sec- 
tion, with arbitrary circulation for any angle of attack 
a. Equation (39) represents the important special 
case in which the circulation is specified by the Kutta 
condition. The various symbols are functions only of 
the coordinates {x, y) of the airfoil boundary and ex- 
pressions for them have already been given. In Tables 
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I and II are ^iven numerical results for different air- 
foils, and explanation is there made of the methods of 
calculation and use of the formulas developed. 

We hav^c immediately hy equation (3) the value of 
the pressure p at any point of the surface in terms of 
the pressure at infinity as 



Some th(‘oretical pressure distribution c\irves are given 
at the end of this report and comparison is there made 
with exi)erimental results. These comparisons, it wdll 
be seen, within a large range of angles of attack, are 
strikingly goo<I.‘^ 

GENERAL WING-SECTION CHARACTERISTICS 

The remainder of this report will be devoted to a 
discussion of the parameters of the airfoil shape affect- 
ing aerodynamic properties wdth a vieW' to determining 
airfoil shapes satisfying preassigned proy)crties. This 
discussion will not only furnish an illuminating sequel I 





Kk.I'HK Ifi. 

to tl\e foregoing analysis leading to a number of new 
results, but will also unify much of the existing theory 
of the airfoil. In the next section we shall obtain 
some expressions for the integrated characteristics of 
the airfoil. We start with the expressions for total 
lift and total moment, first developed by Blasius. 

Blasius’ formulas. — Let C in Figui e 16 represent a 
closed streamline contour in an irrotational fluid field. 
Blasius' formulas give expressions for the total force 
and moment experienced by C in terms of the complex 
velocity potential. They may be obtained in the fol- 
lowing simple juanner. We have for the total forces 
in the x and y directions 

” f pTds= ~ f pdy 
c c 

Pv = y'pi/d.s’ - = J'pdx 
c c 

- iPy = - f p {dy + idx) 
c 

A pa{>er rlevote<l to riiore extensive fxpi>li<‘aUoiis to preseiitMlay airfoil.s is in 
progress. 


The pressure at any point is 


P.-iPv-^fv^idy^-idx) 

^ c 

ip p dtc dw -j- 

= 2id7d?‘*" 

where the bar denotes conjugate complex quantities. 
Since C is a streamline, zvdy - lylz = 0. Adding the 
quantity 

i P f + h\) {ihdy - r ^di-) = 0 
c 


to the last equation, we get. 


Pz -iPv= f {Vj ~ ii\y{dx + idy) 

- c 




The differential of the moment of the resultant 
force about the origin is, 
d A/o = p{x dx + y dy) 

- K. P, of p[x ^ + T/ dy L i (t/dj ~ zdy)] 

^-R. P. oi p z dz 

where P. of” denotes the real part of the complex 
quantity. We have from the previous results 

Then dA/o=-«. P. oi z Az 

and P. ol^^f(^^zdz (43) 

Lot us now for completeness apply these formulas to 
the airfoil A in the f plane (fig. 14) to derive the Kutta- 
Joukow’sky classical formula for the lift force. By 
equation (32) we have 

d^ 2ir2 Z^ 

and by equation (25) 

/I ^ - 2 .,.3 . - • 


d w _ 

di" d^ df 

1'^ ('f Hlasiiis, H: Zs. f. MiUh. n. Phys. Hd. .‘>8 S. *J3 and Bd. 59 S. 43, 1910. 
.‘Similarly, 




a less convenient relation to u.ie than (42). 

Note that when the region about Cis regular the value of the integral (42) remains 
unchanged by integrating about any other curve eTiclosing C. 
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and 



where 



"2 tldr/ dz 


! tion of origin as may he readily verified by employing 
equations (26) and (32') and integrating around the C\ 
circle in the plane. It is indeed a remarkable fact 
that tlie total integrated characteristics, lift and loca- 
I tion of lift, of the airfoil depend on so few parameters 
I of the transformation as to be almost independent of 
i the sha])e of the contour. The parameters R, Oi, 

' and Cl involved in these relations will be dis(‘ussed in a 
later paragra])h. 

I We shall obtain an interesting result by taking 
; moments about the point f = Ci instead of the origin. 
{M ill fig. 17.) By equation (25) we have, 



and by C(] nation (43) 




= - 16 '“ plT 
Therefore 

I\ = pVV sin al 
Py ~ pVV cos a 1 

and are the components of a force pW which is per- 
pendicular to the direction of the stream at infinity. 
Thus the resultant lift force experienced by the airfoil is 

L = pVV (44) 

and writing for the circulation T the value given by 
equation (35) 

L = 4:irRpY^ sin (a + ^) (45) 

The moment of the resultant lift force about the 
origin 0 is obtained as 

M, ll. F. of-|/(gfYf dr 

= 7f. P. + . -)x 

('.+•-+"'4- -x-?" 

= R. P. of — ^2x1 (coefficient of e'^) 

= R. P. of — ^ 27rl {A^ r 2I1C1) 

or, Mo is the imaginar\" part of 7rp(/h + /l,Ci). After ! 
putting Cl ^ me'^ and «i we get 

Mq 2 ttpVW sin 2 (a + 7 ) 4 - pIT m cos (a + 5 ) ( 40 ) 

The results given by equations (44) and (46) have 
physical significance and are invariant to a transform a- 

» It may be recalled that c\A— f and + etj. (2.*/).) 

* n ^ 



This result could have been obtained directly from 
equation (40) by noticing that pVV in the second term 
is the resultant lift force L and that Lm cos (a f6) 
represents a moment which vanishes at I\f for all values 
of a. (In fig. 17 the complex coordinate of M is 
me^\ the arm 01! is m cos (a + 5).) The jierpendic- 
ular JiM from M onto the resultant lift vector is simply 
obtained from 


as 


, sin 2 (a 4 7) 

'^^“27? sin (a+/3) 


(48) 


The intei'section of the resultant lift vector ^\'ith the 
chord or axis of the airfoil locates a point which may 
be considered the center of pressure. The amount of 
travel of the center of pressure with change in angle 
of attack is an important characteristic of airfoils, 
especially for considerations of stability, and will be 
discussed in a later paragraph. 


First obtained by R. von Miscs. (Refereiiee 6.) The work of \ on Mises forms 
an elegant geometrical study of the airfoil. 
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The lift force has been found to be proportional to 
sin {a + /S) or writing? a + = ai 

L = iirpRV^ sin ai (49) 

where may be termed the absolute an^le of attack. 
Similarly writing a +7 = a 2 

Trb^p sin 2 (oO) 

With von Mises (reference 6 , Ft. II) we shall denote 
the axes determined by passing lines through A/ at 
angles jS and 7 to the x axis as the first and second axes 
of the airfoil, respectively. (Fig. 18.) The directions 
of these axes alone are important and these are fixed 
with respect to a given airfoil. Then the lift L is 
proportional to the sine of the angle of attack with 
respect to the first axis and the moment about M to 


If this moment is to be independent of a, 
cients of sin 2a and cos 2a must vanish. 
Then 

6 ^ cos 2y~ Rr cos (jS i- o-) 


and 


sin 2y ^Rr sin(|3 + a) 


the coeffi 


Hence, 



(j^2y— ^ 


Then if we move the reference point of the moment to 

a point F whose radius vector from A/ is the 

moment e.xisting at F is for all angles of attack con- 
stant, and given by 

M 1 ^ = 2 irp¥ sin 2 ( 7 -/?) (51) 



the sine of twice the angle of attack with respect to 
the second axis. 

From equation (47) we note that the moment at any 
point Q w’hose radius vector from M is given by 

A/q = 2TpbW^ sin 2 (a + 7) — Dr cos (a + c) 

Let us determine whether there exist particular 
values of r and u for w^hich Mq is independent of the 
angle of attack a. Writing for L its value given by 
equation (45), 

Mq-2'Kp¥V^ sin 2{a F 7 ) — 4xpi?r sin (a + j3) cos (a + (t) 
And separating this trigonometrically 

A/q- 27 rpF^[( 6 ^ cos 27 — Rr cos (/3+ (t-)) sin 2 a 
+ (b^ sin 27 — sin (/9+<r)) cos 2 a 
— JF?r sin (iS — cr)] 


It has thus been showm that with every airfoil pro- 
file there is associated a point F for which the moment 
is independent of the angle of attack. A change in 
lift force resulting from a change in angle of attack 
distributes itself so that its moment about Fis zero. 

F’rom equation (47) it may be noted that at zero lift 
(l^ e., a = - / 3 ) the airfoil is subject to a moment couple 
which is, in fact, equal to My. This moment is often 
termed “diving moment “ or “moment for zero lift.” 
If My is zero, the resultant lift force must pass through 
F' for all angles of attack and we thus have the state- 
ment that the airfoil has a constant center of pressure, 
if and only if, the moment for zero lift is zero. 

The point F, denoted by von Mises as the focus of 
the airfoil, will be seen to have other interesting prop- 
erties. We note here that its construction is very 
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simple. 1 1 lies at a distance from M on a line makinj^ 

an^le 2y- ^ ^\^th respect to the / axis. From Fifjure 
18 we see that the ang;le }>etween tliis line and tlie first 
axis is bisected by the second axis. 

The arm from F onto the resultant lift vector L 
{hjr is designated FT in Figure 18; note also that 
being perpendicular to Z, must be parallel to the direc- 
tion of flow; the line T\^ is drawn parallel to the first 
axis and therefore angle V 7 'F a - is obtained as 

A _ ^ 

^ L 2/^ sin(a I 
or setting h = sin 2 - 7) 

. _ h 

sin (a I /3) 

But hp is parallel to the direction of a, and the relation 
^ ^ sin (q;+/ 3) states then that the ]>rojection of 
hp onto the line through F perpendicular to the first 
axis is equal to the constant //. {h is designated FV in 
the figure) for all angles of atta(>k. In other words, 
the pedal points T determined bv the intersection of 
hp and L for (ill positions of tln^ lift vector L lie on a 
straight line. (The line is determined by T and V in 
fig. IS.) The parabola is the only curve having the 
property that pedal points of the perpendiculars 
dropped from its focus onto any tangent lie on a 
straight line, that line being the tangent at the vertex. 
This may be shown analytically by noting that the 
equation of L for a coordinate system having F as 
origin and FI' as negative x axis is 

X sin «i -j y cos - hp~ - - 

By differentiating with respect to « f ^ and elim- I 
inating ai we get the equation of the curve which the | 
lines L envelop as y~---Ah{x r-h). From triangle FVS 
in Figure 18, it may be seen that the distance 

A/F— is bisected at S by the line T\ \ for, since 
FI ^ ^ 2 - 7), then 

== 2^* It has thus been shown that the resultant lift 

vectors envelop, in general, a parabola whose focus is 
at F and w hose directrix is the first axis. The second 
axis and its perpendicular at A/, it may be noted, are ^ 
also tangents to the parabola being, by definition, the 

resultant lift vectors ftir o ” 7 and — 7, respec- 

tively. 

If the constant h reducc*s to zero, the lift vectors 
reduce to a pencil of lines through F. Thus a constant 
center of pressure is given by h =-0 or sin 2(/3~7)- 0 
w^hich is equivalent to stating that the first and second 
axes coincide. The lift parabola opens downward 
when the first axis is above the second axis it 

reduces to a pencil of lines when the tw-o axes are I 


coincident (/if ^7) and opens upw^ard w^hen the second 
axis is above the first (/3<^7). 

W. Muller introduced a third axis wdiich has some 
interesting pro])erties. Defining the complex coordi- 
nate as tlu‘ centroid of the circulation by 

and using eejuatiems (25) and (32) one obtains 
where 

si" “ + ^ si" (“ - -^)l] 

^" = 2 “ - 7 ? ''"S (a 4 27)1 1 


The eipiation of the lift vector lines referred to the 
origin at A/ and / axis drawn through A/ is 


X cos a — 7/ sin a = 


F sin (a + 7) 
2 R sin (a+ 


(54) 


and it may be seen that the point (oi, y^) satisfies this 
equation. A'he centroid of the circulation then lies on 
the lift vectors. liy elimination of a from equation 
(53) one finds as the locus of {xq, y,)) 

2 xu{R cos ^- j^ cos (/3-27)] -r27/o[F sin 

(55) 


whicli is the equation of a line, the third axis, and 
proves to be a tangent to the lift parabola. Geomet- 
rically, it is the perpendicular bisector of the line FF' 
joining the focus to the ]>oint of intersection of the 
first axis wuth the circle. (Fig. 18.) 

The conformal centroid of the contour. — It has 
already been seen that the i)oint A/ has special inter- 
esting propei’ties. The transformation from the air- 
foil to the circle having A/ as center was exjiressiul in 
the normal form and permitted of a very small €(s^) 
curve. (See j). 188.) It was also shown tliat the 
moment with res])cct to A/ is simply pn>])ortioTial to 
the sine of twice the angle of attack witli respect to 
the sec(ui(l axis. AVe may note, too, that in the [ires- 
entation of this nqiort the coordinate of A/, F'f'i 



of the rp(ip) curve. 

Wc shall now' obtain a significant pro])crty of A/ 
invariant with respect to the transformation from air- 
foil to circle. We start with the evaluation of the 
integral 


A 'dfi 


« Hefi*rence 7, p, ,\lso Zs. fiir .\un. .Matfi. u. Meoh. IPI, 3 S. 1 17, iy23. 
Airfoils having ih? same Tirst, second, ami third axes are alike llieorelically in 
lotal lift properties and also in travel of the eenier of pres.snre, i. e , they have Ihe 
same lift parahnla. 
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where A is the airfoil contour, the differential of 

arc aloni^ A, and ? as will be recalled, is the 11012 : 111 - 

fication factor of the transformation mapping 

airfoil into circle; i. e., each element ds of A when 

I dz I 

magnified by gives diS the differential of arc in the 
plane of the circle, i. e., \dz\. Then we have, 

I 

/rhiLid.'? - y'r(s) Id^l ami by equation (25), 

A c 

^ z-\- ^ . ."^|<b| 

27T / \ 

= / (p, 

- 2tt R Ci 

Then 


A 

dz] 

<Tf| 

ds 

f 

A 

dz\ 

dr 

|d.9 


The point M of the airfoil is thus the conformal cen- 
troid obtained by giving each clement of the contour 
a weight equal to the magnification of that element, 
which results when the airfoil is transformed into a 
circle, the region at infinity being unaltered. It lies 
within any convex region enclosing the airfoil contour.^'* 


ternal to the line 4a in the f plane maps uniquely into 
the region external to the circle \z\=a. A point Q of 
the line corresponding to a point P at is obtained 
by simply adding the vectors a{eJ ^ or completing 
the parallelogram OPQP' . 


For ^ = 0, we have from equation (6) 
/ = 2a cosh cos 0 = 2a cos 0 
y = 2a sinh ^ sin 0 = 0 


Then the parameters for this case are R^a, /3 = 0, 
a, =a“ (i. e., b = a, 7 = 0), and M is at the origin O. 
Taking the Kutta assumption for determining the 
circulation we have, 


the circulation, r = 47raFsiii a 

the lift, L = 47rapV^^ sin a 

moment about d/, = 27raVI^^ sin 2a 

h '^ . 

position of /'is at Zf = Ci + ^^=a 


(•^7) 


Since j3 = 7 , we know that the travel of the center of 
pressure vanishes and that the center of pressure is at 



ARBITRARY AIRFOILS AND THEIR RELATION TO 
SPECIAL TYPES 

The total lift and moment experienced by the air- 
foil have been seen to depend on but a few parameters 
of the airfoil shape. The resultant lift force is com- 
pletely determined for a particular angle of attack by 
only the radius R and the angle of zero lift /3. The 
moment about the origin depends, in addition, on the 
complex constants c, and or, what is the sqme, on 
the position of the conformal centroid M and the focus 
F. The constants Ci and a, were also shown (see foot-, 
note 20) to depend only on the first and second har- 
monics of the €(<p) curve. Before studying these ^ 
parameters for the case of the arbitrary airfoil, it will 
be instructive to begin wdth special airfoils and treat 
these from the point of view^ of the ‘‘conformal angular 
distortion’' curve. 

Flow about the straight line or flat plate. — As a 
first approximation to the theory of actual airfoils, 
there is the one which considers the airfoil section to 
be a straight line. It has been seen that the line of 
length 4a is obtained by transforming a circle of radius 

a^ 

a, center at the origin, by + region ex- 

P. Frank ami K. I.owiier, Math. Zs. Hd. 3, S. 7S, Al.so reference 

p. H«. i 


F or at one-fourth the length of the line from the lead- 
ing edge. The complex flow potential for this case is 

where — aMs the inverse of equation 

(5). Since \p{(p) = €{<p) ^0 for this case, equation (39) 
giving the velocity at the surface reduces to 


in(| + a) 
r=V\ A2 1 

. <p 

Sin 2 j 
and hy equation (41) v= V 


for r = 4iraF sin a, 

^ pin (v?^q)\ 

\ sin <p 7 


for r = o. 


Flow about the elliptic cylinder. — If equation (5) 
is applied to a circle with center at the origin and 
radius ae'^, the ellipse (fig. 19) 

- - + ....= 1 

(2a cosln/ )“ (2a sinh^)^ 

is obtained in the f plane and the region external to 
this ellipse is mapped uniquely into the region external 
to the circle. The same transformation also trans- 
forms tliis external region into the region internal to 
the inverse circle, radius ae~'^. We note that a point 
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Q of the ellipse corresponding to P at is 

obtained by simply completing the parallelogram 
OPQF^ (fig. 19) where P' now terminates on the circle 
ae~^. The parameters are obtained as R = ae*y /? = 0, 
at = 02 , A/ is at the origin 0. Then, assuming the rear 
stagnation point at the end of the major axis, 

r = 47rae’^F sin a 
L = 4TTf>ae'f'V^ sin a 
A/jif = 27raVV"^ sin 2a 


Since /3 = 7 , the point Pis the center of pressure for all 
angles of attack and is located at from O or a 

distance ae* from the leading edge. The quantity 


EF ae'^ 
EE'~2a{e'f' + e-'i') 


cosh \p + sinh \J/ 
4 cosh ip 


^ (1 +tanh 


represents the ratio of the distance of F from the 
leading edge to the major diameter of the ellipse. 

The complex flow potential is identical with that 
given by equation (58) for the fiat plate, except that 
the quantity in the numerator of the second term is 
replaced by the constant Since \p{(p) = constant, 

€(v?) = 0 and equation (39) giving the velocity at each 
point of the surface for a stagnation point at end of 
major axis becomes 

.._T >in (y+g)+sin a]e^ ' 
and for zero circulation by equation (41) ' 

( 59 ') 

Vsinh ^^ + sin^^ 


Circular arc sections. — It has been shown that 

the transformation applied to a circle with 

center at 2 = 0 and radius a gives a straight line in the 
^ plane, and when applied to a circle with center 2 = 0 
and radius different from a gives an ellipse in the f 
plane. We now show" that if it is used to transform a 
circle wi th cen ter at z = is {s being a real number) and 
radius + a circular arc results. The coordinates 
of the transform of the circle C in the ^ plane are given 
by equation (6) as 

X = 2a cosh \j/ cos 6 
y = 2a sinh ^ sin 6 


A relation between \j/ and 6 can be readily obtained. 
In right triangle OMD (fig. 20), 0M=s, angle OMD = d, 
and recalling that the product of segments of any 
chord through O i.s- equal to a\ 01)=}^ (OP~OPi)^ 




= a sinh yp. 


Then s sin 0 = a sinh and from 


the equation for y,y = 2s sin^0. Eliminating both 6 and 
\p in equation (6) we get 




(60) 


the equation of a circle; but since y can have only 
positive values, we are limited to a circular arc. In 
fact, as the point P in Figure 20 moves from A' to A 
on the circle, the point Q traverses the arc Ai' Ai and 
as P completes the circuit A A' the arc is traversed in 
the opposite direction. As in the previous cases, we 
note that the point Q corresponding to either P or to 
the inverse and reflected point P' is obtained by com- 
pleting the parallelogram OPQP\ We may also note 


y ■ / Curve 



Figure 20. — The circular arc airfoi] 

that had the arc AiAi' been preassigned with the 
requirement of transforming it into the circle, the most 
convenient choice of origin of coordinates w^ould be 
the midpoint of the line, length 4a, joining the end 
points. The curve B then resulting from using trans- 
formation (5) would be a circle in the 2 ' plane, center 
at 2 ' = is, and the theory developed in the report could 
be directly applied to this continuous closed B curve. 



Had another axis and origin been chosen, e. g., as in 
Figure 21, the B curve resulting would have finite 
discontinuities at A and A\ although the arc AjA/ is 
still obtained by completing the parallelogram 0PQP\ 
The parameters of the arc AjA/ of chord length 4a, 
and maximum height 2s are then, 

^ = tan~^^* The focus F may be constructed by 
erecting a perpendicular to the chord at A' of length s 
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and projecting its extreinity on .1/^1'. The center Al' 
of the arc also lies on this line. 

The infinite sheet having the circular arc as cross 
section contains as a special case the flat plate, and 
thus permits of a better approximation to the mean 
camber line of actual airfoils. The complex flow poten- 
tial and the formulas for the velocity at the surface 
for the circular arc are of the same form as those 
given in the next section for the Joukowsky airfoil, 
where also a simple geometric interpretation of the 
parameters € and \[/ are given. 

Joukowsky airfoils. — If equation (5) is applied to 
a circle with center at s being a real number, and 
with radius R — a symmetrical Joukowsky air- 

foil (or strut form) is obtained. The general Joukow- 
sky airfoil is obtained when the transformation 

^ = 2 + ^ is applied to a circle C passing through the 

])oint z~—a and containing z-^a (near the circum- 
ference usually), and whose center A/ is not limited to 
cither the xory axes, but may be on a line OM inclined 
to the axes. (Fig. 22 .) The parametric equations of 
the shape are as before 

X = 2 a cosh yp cos d 

. . . (6) 

y = 2a sinh p sin 6 

Geometrically a point Q of the airfoil is obtained by 
adding the vectors and or by completing 

the parallelogram OPQP' as before, but now P' lies on 
another circle B' defined as the inverse 

and reflected circle of B with respect to the circle of 
radius a at the origin (obtained by the transformation 
of reciprocal radii and subsequent reflection in the x 
axis). Thus OP* OP' = a* for all positions of P, and 
OP' is readily constructed. The center A/i of the 
circle B' may be located on the line AM by drawing 
OA/j symmetrically to OA/ with respect to the y axis. 
Let the coor<linatc of A/ be z ^ is ^ where d, s, 
and d are real quantities. The circle of radius a, with 
center A/n at 2 = is, is transformed into a circul/ir arc 
through A^A^' which may be considered the mean 
camber line of the airfoil. At the tail the Joukowsky 
airfoil has a cusp and the upper and lower surfaces 
include a zero angle. The lift parameters are 

R~ + (/, tan~* ai = a^ = or 6 = a and 

7 = 0. Since 7 = 0, the second axis has the direction of 
the X axis. The focus F is determined by laying off 
a} 

the segment AfP=^ on the line A//1'. This quantity, 

it may be noted, is obtained easily by the following 
construction. In triangle A1DC\ MD^R, AIC' and 
MC are made equal to a, then CF drawn parallel to 

DC' determines MF=-^’ The lift parabola may be 

now determined uniquely since its directrix AM and 
focus F are known. 


It may be observed that if it is desired to transform 
a preassigned Joukow-sky profile into a circle, there 
exists a choice of axis and origin for the airfoil such 
that the inverse of transformation (5) will map the 
airfoil directly into a circle. This axis is very approx- 
imately given by designating the tail as ( — 2 a, 0 ) and 
the point midw’ay betw^een the leading edge and the 
center of curvature of the leading edge as (-r 2 a, 0 ) the 
origin then bisecting the line joining these points. 



Fi<.l KE 22.— The Joukowsky airfoil 


The complex potential flow^ function for the Jou- 
kowsky airfoil is 

w(f) = - + ;;({■) (01) 

w^here 

By equation (39) we have for the velocity at the 
surface 

yisin(a+ if) -T sin(a+ 

r_.__ 

y (sinh V + sin=9)( 1 + ( ^) ) 

This formula w as obtained by transforming the flow^ 
around C into that around B and then into that around 
A. Since w'e know that B is itself a circle for this 
case, w^e can simply use the latter two transformations 
alone. 

We get 

V[sm + if) (a+ 

Vsinh ^^ + sin^0 

That these formulas are equivalent is immediateh^ 
evident since the quantity 



is unity being the ratio of the magnification of each 
arc element of C to that of B. (See eq. (37).) 
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A voiy siiii])Io ti ichI piotioo of tl:o |:an n elors 

t mul \l/, exists for the east’s discussed. In Figure 
23 the value of t or at the i)oint I* is siiii])!y 


1 

The constant - J'\pdip is readily shown to be 

27T 



angle 0PM, i. e., the angle sii]> tended at P t)v the 
origin O and the center M. The angle of zero lift is 
the value of e for i e, cyan "" i3 "" d/. In 

particular, we may note that e 0 at S and S\ which 
are on the straight line OM. Consider the triangle 

OMl\ where OP^ae*, MP R--ae*% aii-le 

OPM=t; also, MOX^S, MOP -6-5, OMP -^r 
(v?“5). Then by the law of cosines, we have 

\ ... Op ^ ] f p-’ 

or 

\ h)g (1 4 2p cos(<^ — 5) -r p^) (<)3) 




and bv the law of sines 


p sin((p 8 ) 

( I -r- 2p cos( <p (5) ■ p'^)''^ 


, . p Sin(<p 6 

e(«pr tan 

1 p (*OS((p d) 


(b4) 


We see that, as re(|uired, the expressions for the 'T'adial 
distortion’' i/'(ip) and the '‘angular distortion’’ e(<p) 
are conjugate Fourier series and may be expressed as 
a single complex fpmntitv 

“■ ^ 

I 

=-log [L -f pe’*^'""^>] 

It is evident also tliat the coeilicient for = 1 or tlie 
‘‘first harmonic term” is simply and a translation 
by this quantity brings the (‘iicle C into coincidence 
with B as was pointed out on i)age 1S7. 


invariant to the clioicc of origin 0, as long as 0 is 
within B. We have 


j 2tT ^ 27T j 

f^pd^= f . log (I -T 2p cos (<p-5j -p'je"Vi<P 
“^0 “''U “ 



00 


<A.. i i; (-1) 

1 


n - 1 


C()S 


n(ip-6) 

n 



dip = \pii 



^{ip), i(d), and e(p') curv'es for this airfoil. 



Arbitrary sections. — In order to obtain the lift 
parameters of an arbitrary airfoil, a convenient choice 
of coordinate ..aCs is first made as indicated lor the 
Joukowsky airfoil and as stated previously. (Page 181.) 
The curve resulting from the use of transformation (5) 
will yield an sirhitrary curve which will, in 

general, differ very little from a circle. The inverse 
and reflected curve will also be almost circular. 

The transition from the curve to a circle is 

reached hy obtaining the solution c-(<p) of equation 
(13). The method of obtaining this solution as 
already given converges with extreme ra])idity for 
i nearly circular curves. 
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The picture is analoj^oiis to tlmt ^iv'cti 

for Uh‘ s])C(‘ial ruses. Jri Fi^nire 20 it inuy he seen that 
a point Q on the jiirfoil (X. A. C. A. -MOJ corre- 



sponding to I* on th(‘ B ('urv(‘ (or P' on tlu‘ fi' (Mirve) 
is obtained by constructing^ ])a.rallelo^raiu OPQP' . 
The yl/{0) and eiO) cui‘v(‘s are shown in Fijj^ui'e 27 lor 
this airfoil. I'he (‘onipl(‘x velocity ]>otential and the 
expn'ssion foi’ velocity at the surface arc j^iven resjx'c- 


The method used for arbitrary airfoils is readily 
a])j)lied to arbitrary thin arcs or to broken lines such 
as the sections of tail surfaces foi’iu approximately. In 
Kiirnre 2() the ])art of tin' aii'foil boundary ab(»V(‘ the x 
axis transforms by c(iuation (o) int<» the two discon- 
tinuous arcs sliown by full lines, while the lower 
boundary transforms into the arcs shown by dashed 
lines. If the u])])er boundary surface is alone ^iven 
(thin airfoil) we may obtain a (“losed curve only 

by joining tlu' end ])oints by a chord of leiiirth 4ri and 
choosing the ori^dn at its midi)oint.‘^ The resulting 
(‘urv(‘ lias two double ])oints tor whiidi the lirst deriva- 
tive is not uni(piely <letined and, in general, it may be 
seen that infinite velociti(‘S correspond to smdi ])oints. 

At a ])oint. of the 4^{6) (uirve <‘(»rrespondin^ to a 
matluMuatically sharj) coi ruu-, there exist two taiu^ents, 

that is, the sloiie is finitely discontinuous. 'Fhe 

' ^ A a 



tively l)y e(juations (33) and (30). The lift param- 
eters are 

R 

R = ^ ^ fzvni (at 0 7r),A/is at£: - Ci - J' 

^ 0 

and /As at :: Cj ^ where a, is pv(m in e<iuation (25^). 

The first and second axes for the XX A. A. -Mb 
airfoil are found to ci)incide and this airfoil has then a 
constant ccmter of pressure at F. Fij>:ures 2<S (a) to 
28 (1) lh(‘ )uessure distribution (alon^ the / axis) 
for a series of an^h‘s of attack as calculated by this 
tiu'orv and as obtained by ex])eriment."^ Table 1 
contains the essential numerical data for this airfoil. 

** 'I’he ex[)oriitienlal rf.MiUs nri? taken from test No. 32.'l of tlie N', A C, variable- 
(iensity wital tnntiel. The angle of attack » substituted in eriiiation (W) ha.s lieen 
inoflificd arbi(raril.v to lake account of the ellecls of finite simn, tunnel-wall inter- 
ference, and visco-sity, by choosiriR it so that the theoretical lift is aljont 10 per cent 
tiKjre than the (.(jrresporiding experimental value. The actual values of the lift 
( rjelUcient.s are given in the figure.s. 

4070S ::4-- 14 


: curve i{9} must liavo an infinite slope at siicli a ])oint 
I for according to a tlieorem in the theory of Fourier 
I scries, at a ])oint of diseoiitinuity of a F. S., the eon- 
I jupite F. S. is pro])erly divoiTcent. This manifests 
I itself in the velotuty-fonnula equation (39) in the fac- 

I tor ^ 1 f which is infinite at these sharp corners. 

For ])raetical purposes, however, a roundinjx of the 
! sliarp edge, however small, considerably alters the slope 

I at this ])oint. 

I do ^ 

I Ideal angle of attack.~~A thin airfoil, represented 
■ by a line arc, has both a shar]> loading edge and a 
shar]) trailing edge. The Kutta assumption for fixing 
^ the circulation places a stagnation point at the tail for 
' all angles of attack. At the leading edge, however, 

Note that 7 («+jt> = -T(0) for this case. 
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Jheoreftca! 



Experimenfal 

o Upper surface x Lov^er surface 

(A verage R.N. = 5 x tO^} 






Fioi'RES 2S a to e.— Tlieorelical and exjwriniental jiressure distriliutioa fur tlie M 6 airfoil at various angles ol attack 
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tlio velocity is infinite at all aiiirk's of attack except 
one, namely, that an^le for which tlie other stagnation 
])oint is at tlie leading edjzo. It is natural to expect 
that for this anisic of attack in actual cases the fric- 
tional losses are at (»r Tunir a niiniimitn and Urns arises 
the conce])! of ‘‘ideaP' an^k' of attack introduced hy 
Theodorscu (reference 8) and which lias also been 
designated “anjrle of best sl r(\ainlininp:/’ Tlie defini- 
tion for the ideal anjrle may be extended to thick 
airfoils, as that anirle foi‘ which a sta«;nation jioint 
occurs directly at the foiimiost ])oint of the mean 
camber line. 

The lift at the leadini: vaiiislies and the chan<;e 

fiom velocity to pressure al(ui^ the airfoil surface is 
usually more gradual than at any other angle of attack. 


of this functi^m, one can determine airfoil sha]>es of 
(hdinite ])n^])(‘rti<‘s. The €(<^) function, which we have 
designated conformal angular distortion function, will 
be seem to diUermine not only the shape but also to 
give easily all the tlu'oretical aerodynamic charactcu'- 
isti(‘s of the airfoil. 

An arliitraiy e(v) (uirve is cho.sen, single valued, of 
]>ei'iod 2tt, of /(M‘o area, and such that -- oo < 


These limiting 


values of ai’c 
dv? 


r ^ 1. 


far bevond value: 


Theory fiC3l 


Experimen fa/ 

O Opper surface 
X /.oiver surface 


yielding airfoil shapes/'" The function, excejit for 
the constant is given by the conjugate of the 
Fourier expansion of €iV) or, what is the same, by 
evaluating eij nation (14) is a deliiiite integral. The 

" •CG'/ Experimer fa! 

O Upper surface 
X Lower surface 


~4.0 

-3.5 

-3.0 

-35 

-30 


-(.5 


£ 0 
7 


5 

J.O 



The minimum profile drag of airfoils actually occurs 
very close to this angle. At the ideal angle, which we 
denote by a/, the factor [sin (a I v?) I sin (cy [ /3)] in 
eijuation (39) is zero not only for 0 - tt ore tr " ^ but 
also for d ---- 0 or e €,v. We get 

(Xj I 6,v " : 6'/’) or 


CREATION OF FAMILIES OF WING SECTIONS 

The ])rocess of transforming a cind(‘ into an airfoil is 
inherently less <Iillicult than the inverse process of 
transforming an airfoil into a circle. By a liirect ai)]>li- 
cation of previous n‘sults we can dei’ive a powerful and 
flexible method for the creation of general families of 
airfoils. Instead of assuming that the i^(^) curve is 
])reassigned (that is, instead of a given airfoil), we 
assume an arhitrary or €()/) curve as given. 

This is eipii valent to assuming as known a houndary- 
value function along a lun le and, by the {iroper choice 

SiiltjtH't to some reslrio iors L'iveti ii; ttio iieU jiarui'rapti. 


constant is an important arbitrary parameter 
whi(*h ])crmits of (dianges in the shape and for a certain 
range of vtihu's may determine tJie sharpness of tlie 
trjuling t'dg(*. 


I We first obttiin tlu‘ variable 9 as 9 (ip) «/ — f (sp), i^o 
; tliat the quantity 4^ considered as a function of 9 is 
! 4/ {9) 4^[ip (6/)], The coordinates of th(‘ airfoil surface 

an* then 

X 2<i cosh 4^ cos 9 


y 2a siidi 4^ sin 6. 


1 >) ) 


For coiniiio': lurfuils, w itli a prtip.'r rhoire of orij:in, 

For romiiio:! aitf.iils 4-t is usually hiaweeu 0 <)*j an<l fU*. 'Hu.* ronstant is 
jiut. howysiT, ( urhUniry. Wo havt; li.n*!) tliui the rorulitiori driven liy 
ip2; i:i auHii ieKt to yiel'i ii t'ontrrir free from douttle |>uiut> in the 2 ' pkine. 
' We may also <tu(e r!ie criterion that the itiverse of P(|uatio:i (.^j aioilieh to this 
I eontoiir shall \ i.'irl a roalour in the C plane free from iloiihk* points. ('luisiihT the 
ftiiu'tion for d vary i:ttt from 0 (0 r o ily. The iii*L'ative of earh value nf ^{6) i.i 
this rant-'e i.s rotisidereU as.soiiate'i with —9. i. e., Deseitralc the fiiticlion 

ilius formed from 9--0 to 2ir by Then represents a line are in the r 

I plane, i. e., tlie upper surface of a eonioiir. iSee footnote 'J.',] Then for Hie entire 
I I'fuUour to be fn‘i' from double ixhnts it i.s ne.-e.ssary Unit the lower surface .should iiol 
cross the u]>per, tool i.-. The or'ijincu curve for 0 varying' from w to 2r must not 
cross below 
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'•'lie velocity at the sJirface is 

.. V^_., [sin (g + <p ) 4 - sin ^ 

V "'f' [( ^ ) + ( il t ) ] 

and is obtained by usin^ equation (370 instead of (37) 
in d(‘rivin^ (39). The an^le of zero lift jS is ^iven by 
V? {0) 0 • € {$) for d ■- 7T, i. e., (tt) = tt /3. 

Tbe following figures and examples will make the 
proc(‘ss clear. We may first note that the nu)st natiiral 
method of s])ecifyintt tlie e (v?) function is by a Fourier 
series expansion. In this sense then the elementary 
types of t (v?) functions are the individual terms of 
this ex])ansion. 


I 29(h) to 29(t). In particular, the second harmonic 
term may yield shapes, and by a proper combina- 
tion of first and second harmonic terms, i. e., by a 
proper choice of the constants /li, A2, 61, and 62 in the 
relation 

€(y)=/li sin {<p — bi)+A 2 sin (2<,i? — 62) 

it is ])ossiblc to fix the focus F of the lift parabola as 
the center of pressure for all anj?les of attack.*® The 
e< [nation 

€((^) = 0.1 sin (cp — 60°) + 0.05 cos 2^> 
represents such an airfoil and is shown in Figure 29(u). 

The general process will yield infinite varieties of 
contours by superposition of sine functions; in fact, if 


f(ipj=0J sin<f 
ipo-O.i 



\^'o=0 t 



(!<f; -0.05 

Ipr.’^O.O'D 





( (ipj ^0.! S!n(3<p- 90 7 
-O. / 



t(<p)^O l s>-nf^-45y 




€(ip)= 0.075 S!n(4(p) 



f(ip)^0.! sm(ip -GOy 



S:nO<p 



f (^J =0.075 sm(4^-45 7 



€>ipj^-o.i' smf^- 75 y 



s,T)f^-90y 




Consider first the effect of the first harmonic term 
€ (if) =Ai sin (</5-5 i), 

In Figures 29(a) to 29(g) may be seen the shapes 
resulting by displacing 5, successively by intervals of 
15° and keeping the constants ^h^O.lO and i/'o'^^O.lO. 
The first harmonic term is of chief influence in deter- 
mining the airfoil shape. The case f(y) = 0.1 sin 
(<^--45°) is given detailed in Table II. (This airfoil 
is remarkably similar to the commonly used Clark Y 
airfoil.) The entire calculations are characterized by 
their simplicity and, as may be noted, are completely 
free from tfie necessity of any graphical evaluations or 
constructions. 

Tlie effect of the second and higher harmonics as 
well as the constant \po may be observed in Figures 


the process is thought of as a boundary-value problem 
of the circle, it is seen that it is sufficiently general to 
yield every closed curve for which Riemann’s theorem 
applies. 


Lan(;ley Memorial Akhoxautical Laboratory, 
National Advisory Committee for Aeronattics, 
Langley Field, Va., November 4 , 1032. 

” Thi isacfomplished as follows: Weseek to determine tluu-otisljuits . li. i.Ai.and 

r,? 

^2 SO that S-y, where y is obtained from ef|iiation (2.7) as = = bri and 

we inav note that and Aje*h. These relations are transcernien* 

tal; however, with but a few practioe trials, solutions ran be obtaitie<l at will. AthJi- 
tion of higher harmonics will \ iel<l further shapes having the same center of jiressure 
properties if/J is kept unchanged. 
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APPENDIX 


I, EVALUATION OF THE INTEGRAL. 

2w Q Z 


TT Q a<p 


(p"<p 


&<p 


(13) 


(13') 


The function ^(v>) is of period 2tt and is considered 
known. (Note that the variables <p and tp' are re- 
placed by B and B' , pi and pn and <p 2 'j etc., in 
equation (21) and that the following formula is 
applicable for all these cases.) 

A 20-point method for evaluating equation (13) as 


where the interval p' — s to s is small. And, in 
general, 

(p-2 / 

i f 'PM cot A,p 

Vi 

is very nearly 


■ log 


<p ~ ip.^\ 




1 Reference 2, p. 11, gives a 10-point method result. 


where the range <P 2 ~-*Pi is small and 4^a is the average 
value of \p{<p) in this range. The constants a„ for the 
20 divisions chosen above are actually 


sin 7T - 


= log 


2n+l 


40 


SmTT- 


2n-l 


40 


(n- “9, . 


+ 9) 


As an example of the calculation of e(0) we may refer 
to Table I and Figures 26 and 27 for the N. A. C. A. 
— MG airfoil. From the \f/(d) cur\^e (fig. 27) we obtain 

the 20 values of and ^ for 20 equal intervals of 6. 


a definite integral gives 

For the airfoil (fig. 2G) 

we get the following v 

allies: 

£(<p')^-^j^ao— -rai(lAi- iA-i) +«2(\A2~ ^- 2 ) 

(Upper 
6 surface) 



(Lower 
6 surface) 


fl7 

Td 

'i' (’/'j “* '/'-9)^v=<p/ (I) 

0 (nose) 

0. 192 

0. 000 

^ 0. 049 

- 0 . 002 

where 

TT 

To 

. 185 

. 027 

127T 

10 

057 

. 050 

^1 = value of ^(v?) at <,7 = ^'+ 

27T 

lb 

- 192 

. 000 

137T 

10 

071 

. 030 

71 7T 

= value of ^{<p) at <p ---p 4-^^ 

Stt 

To 

. 189 

-. 030 

Htt 

To 

077 

. oil 

(n=l,-l, 2,-2 9, -9). 

4 7T 

To 

. 174 

0G4 

IStt 

10 

079 

. 000 

and the constants arc as follows: ^ = 0.3142; 

5w 

To 

. 14() 

095 

IOtt 

To 

082 

. 016 

a, = 1.091; ao = 0.494; <73 =-0.313; 04 = 0.217; 05 = 0.158; 
06 = 0.115; 07 = 0.0884; Os = 0.0511; and Og = 0.0251. 

Ott 

To 

. no 

114 

1 77T 

nr 

090 

. 039 

This formula may be derived directly from the 
definition of the definite integral. The 20 intervals^ 

7 7T 

TO 

. 077 

080 

IHr 

”10 

111 

. 091 

1 7T , , 7T . 37T 

chosen are ~ to v? r + to etc. 

Stt 

To 

. 052 

-. 066 

1 Ott 
10 

150 

. 154 

It is only necessary to note that by expanding \p{<p) in 

9r 

10 

. 041 

. 025 

27t ( nose) 

192 

. 000 

a Taylor series around p we get 

7T (tail) 

. 055 

. 000 




^ 'PM cot ^-:^d^~-2. 

p'-s 2 L = 

The V 

alue of 

e at the tail (i. e., the angle 

! of zero lift) 

is, for exam})le, 

using formula I 




+ 1.091(.049-.041) 

+ .494C057-.052) 

+ .313(.071-.077) 

+ .217(.077-.110) 

+ .158(.079-.146) 
-f.ll5(.082-.174) 

+ .0884 (.090 -.189) 

+ .0511U11--192) 

+ .0251 (.150-. 185)1 = .0105 
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The value of e for for example, is obtained by 

a cyclic rearrangement. Thus, 

+ 1.091(.174-.192) 

+ .494(.146-.185) 

+ .313(.110-.192) 

+ .217(.077-.150) 

+ .158(.052-.lll) 

+ .115(.041-.090) 

+ .0884 (.055 -.082) 

+ .0511C049-.079) 

+ .0251 (.057 - .077)] = .0347 

The 20 values obtained in this way form the «i(0) 
curve, which for all practical purposes for the airfoil 
considered, is actually identical with the final i{d) 
curve. 


II. NOTES ON THE TRANSFORMATION. 

S ^/(z) “ c, + z+ -' + -)+ . . . (4') 

£1- 


Thore exist a number of theorems giving general 
limiting vahies for the coefficients of the transforma- 
tion equation (4), which are interesting and to some 
extent useful. If ^=f{ 2 ) transforms the external 
region of the circle C of radius R in the 2 plane, into 
the external region of a contour A in the ^ plane in a 
one-to-one conformal manner and the origin ^ = 0 lies 
within the contour A (and /'( ) = 1) then the area S 
inclosed by A is given by the Faber-Bieberbach 
theorem as ^ 




CO 

V 

n 1 




Since all members of the above series term are positive, 
it is observed that the area of C is greater than that 
inclosed by any contour A in the f plane (or, at most, 
equal to the area inclosed by A if is a circle). 

This theorem leads to the following results 

(a) 

\c^\S2R (b) 

Let us designate the circle of radius R about the 
conformal centroid M as center as t7i (i. e., the center 
is at r circle has been called the “Gruhd- 

kreis’^ or “ basic circle by von Mises). Then since 

d I 

^ represents the distance of the focus F from M, the 

relation (a) states that the focus is always within Ci. 
In fact, a further extension shows that if Tq is the radius 
of the largest circle that can be inclosed within A, then 

F is removed from Ci by at least 


* For details of this and following statements see reference 5, p. 100 and p. H7, and 
also reference 6, Part II. 


From relation (b) may be derived the statement that 
if any circle within A is concentrically doubled in radius 
it is contained entirely within a circle about M as 
center of radius 2R. Also, if we designate by c the 
largest diameter of A (this is usually the ‘^chord” of 
the airfoil) then the following limits can be derived: 


4 




These inequalities lead to interesting limits for the 
lift coefficient. Writing the lift coefficient as 




L 


where by equation (45) the lift force is given by 
L = 4iTrRpV^ sin (a+|3) 

we have 

2tt sin (a + iS.) ^ Cl = sin (a + /3) ^ 47t sin (a + /3) (II) 

The flat plate is the only case where the lower 
limit is reached, while the upper limit is attained for 
the circular cylinder only. We may observe that a 
curved thin plate has a lift coefficient which exceeds 
2ir sin (a + /3) b}^ a very small amount. In general, the 
thickness has a much greater effect on the value of 
the lift coefficient than the camber. For common 
airfoils the lift coefficient is but slightly greater than 
the lower limit and is approximately l.lX27r sin 

(ot + ^). 

Another theorem, similar to the Faber-Bieberbach 
area theorem, states that if the equation ^=f{z) trans- 
forms the internal region of a circle in the 2 plane into 
the internal region of a contour B in the ^ plane in a 
one-to-one conformal manner and/^(0) 1 (the origins 

are within the contours) then the area of the circle is 
less than that contained by any contour /i. This 
theorem, extended by Bieberbach, has been used in an 
attempt to solve the arbitrary airfoil,^ The process 
used is one in which the area theorem is a criterion as 
to the direction in which the convergence proceeds. 
Although theoretically sound, the process is, when 
applied, extremely laborious and very slowly con- 
vergent. It can not be said to have yielded as yet 
really satisfactory results. 


in. LOCATION OF THE CENTER OF PRESSURE FOR AN 
ARBITRARY AIRFOIL 

It is of some interest to know the exact location of 
the center of pressure on the x axis as a function of the 
angle of attack. In Figure 30, (7 is the origin, M the 
conformal centroid, L the line of action of the lift 
force for angle of attack a. Let us designate the 

» Mailer, W., Zs. f. angew. Math. u. Mech. Bd. 5 S. 397, I92.‘i. 

Hdhnriorf, F., Zs. f. angew. Math. u. Mech. Bd. 6 S. 265, 1926. 

Also reference 5, p. 185. 
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intersection of L with the x axis of the airfoil as the 
center of pressure P. 


tained graphic all}^ of the ^ against 6, and « against 
6 curves, respectively; (12) is given by 


In the right AONM we have, 

OM = Cl = = Ai + iBi 

ON = m cos 6 = Ai 
m sin 8^Bi 

and in right AJKM, 

’ sm a sin a 

Then 

sm a 

and A^P^KN tun a sec a — Bi tan a ' 

By equation (48) 

h (« + t ) 

L 2i? sin (a 4-/3) i 

Then the distance from tlio origin to the center of j 
pressure P is 

OP = ON + A^P Ai - Bi tan a ' 

2 (g 4 - 7) i 

^2i?cos a sin (a-f/3; (HI) | 



EXPLANATION OF THE TABLES 

Table I gives the essential data for the transform a- , 
tion of the N. A. C. A. -MG airfoil (shown in fig. 26) 1 
into a circle, and yields readily the complete theoretical 
aerodynamical characteristics. Columns (1) and (2) 
define the airfoil surface in per cent chord; (3) and (4) 
are the coordinates after choosing a convenient origin 
(p. 181); (5) and (6) are obtained from equations (7) 
and (8) of the report; (9) is the evaluation of equation 
(13) (see Appendix); (10) and (11) are the slopes, oh- 


(■ 4 ;) 




-y/ (sinh V + sin^ff)( 1 + ) 


where ^o=- 


1 

2tt 


J" ^ {<p) A<p and may be obtained graphi- 
0 


cally or numerically ; column *(13) gives <^ = 0 4-€. The 
velocity for any angle of attack, is by equation (39) 
V ^ Vk [sin (a 4 - V?) 4- sin (a + /3)] 
and the pressure is given by equation (3). The angle 
of zero lift /3 is the value of c at the tail; i. e., the value 
of e for 0= 7T. 

Table II gives numerical data for the inverse process 
to that given in Table I; viz, the transformation of a 
circle into an airfoil. (See fig. 29.) The function 
€(<^)=0.1 sin {(p — 4o^) and constant ^o = 0*10 are 
chosen for this case. Then ^(v?)=0.1 cos (<^-45°) 
fO.lO. It may be observed that columns (11) and 

(12) giving the coordinates of the airfoil surface are 
obtained from equations (6) of the report. Column 

(13) is given b}^ 


\/("i“h [(j - + (0] 

and the velocity at the surface is by equation (39') 


V = Vk [sin (a 4- 4- sin (a 4- /3) ] 
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TABLE I 


N. A. C. A — M6 


uppp:k surface 


Per cent 
c 

j/ in i>er 
cent c 

X 

V 

sin ^0 

sinA V 

0 

radians 

4' 


(te 

d* 

do 

k 


0 

1,25 

2.50 

5.0 

7.5 

10 

15 

20 

30 

40 

50 

60 

70 

80 

90 

95 

100 

0 

1.97 
2. 81 

4.03 
4.94 
5.71 
6.82 
7. 55 
8.22 

8.05 
7.26 

6.03 
4.58 

3.06 
1.55 

.88 

.26 

2 . 037 
1.9H6 
1.936 
1.835 
1.734 
1.633 
1.431 ; 

1.229 
.825 
.421 
.017 
-.387 
791 
-1. 195 
-1..599 
-l.HOl 
-2.003 

0.000 
.0796 
. ii;i5 
. if.:i 
.200 
.231 
.276 
.305 
.332 
. 325 
.293 
. 244 
. 185 
. 124 
.0t)3 
.036 
.000 

0,000 
. 0465 
.0941 
. 187 
.275 
. 357 
.507 
.6.36 
.8.35 
.957 
1.000 
.963 
.845 
.645 
. 36.3 
. 191 
.000 

0. 0373 
.0341 
. 0342 
.0354 
. air.3 
. 037.3 
.0375 
.0366 
.0330 
. 0276 
. 0215 
.0154 
.0101 
.0059 
.0027 
.0016 
.oa30 

0.000 
.217 
.312 
.447 
. 551 
. 640 
.792 
.923 
1. 153 
1.361 
1.571 
1.764 
1.975 
2.209 
2. 495 
2.690 
3. 142 

0. 192 
. 184 
.184 
.187 
. 189 
. 192 
. 193 
. 190 
. 181 
. 165 
. 146 
.124 
.100 
.077 
.0.52 
.040 
.055 

Milo 

0.000 
-.010 
.009 
.022 
. 022 
.020 
-.009 
-.031 
~.0()6 
-.085 
100 

-. 100 
100 
-.088 
-.067 
-.o:i5 
.000 

0. 085 
.080 
.080 
.080 
.080 
.085 
.095 
.100 
. 107 
,088 
.060 
.025 
-.029 
-.056 
-.085 
-.080 
-.027 

6. 280 
4. 249 
3. 368 
2. 557 
2. 163 
1.929 
1.660 
1.498 
1.324 
1.220 
1. 166 
1.152 
1. 167 
1.302 
1.687 
^ 2.340 
\9,S3 

-2 37 
10 52 
16 41 
25 4 

31 31 
37 3 

46 27 
54 38 
69 11 
82 7 

94 55 
106 25 
118 28 
131 19 
146 30 
1.56 28 
180 36 


LOWER SURFACE 


0 

1.25 
2. 50 
5.0 
7,5 
10 
15 


40 

50 


95 

100 


2.017 
l.OHfi 
i.9:«i 
1.835 
1.734 
1.633 
1.431 
1.229 
. H25 
.421 
.017 
-.387 
791 
-1. 195 
-1.599 
-1.801 
-2.003 


0 . 000 
-.071 
I -.089 
-.110 
I -.122 

-. 131 
-. 140 
-. 146 

153 
I -.158 

-.159 

154 
-. 141 
-.114 
-.072 
-.044 

.000 


0.000 
. 0425 
.0844 
. 173 
. 259 
. 342 
.494 
.626 
.831 
.956 
1.000 
.963 
.845 
. 645 
. 363 
. 191 
.000 


I. 0373 
.0297 
.0234 
.0176 
.0144 
.0125 
. (K199 
. 0085 
.0070 
. 0065 
.0063 
. OOf52 
. 0058 
.0050 
.0035 
.0025 
.0030 


6.283 
6. 075 
5. 989 
5. 855 
5. 749 
5, 659 
5.505 
.5. 371 
5. 136 
4, 924 
4.712 
4.518 
4. 307 
4.074 
3. 788 
3.594 
3. 142 


0. 192 
. 172 
. 152 
. 132 
. 120 
. 112 
.099 
.092 
.084 
.081 
.079 
. 0785 
.076 
.071 
.0.59 
.050 
. 055 


-0- 0457 
-.0781 
-.0850 
-.0882 
-.0850 
-.0811 
-.0723 
-.06,37 
-.0516 
-.0421 
-.0350 
-.0310 
-.0300 
-.0295 
-. 0^35 
-.0140 
. 0105 


0. 000 

0.085 

6.280 

-2 

37 

. K13 

. 120 

4.615 

-16 

21 

. 160 

. 050 

3.525 

-21 

43 

. i;i3 

-.010 

2.510 

-29 

35 

, 109 

-.045 

2. 025 

-35 

28 

.080 

-.057 

1.764 

-40 

24 

.069 

'-.067 

1 . 4«J 

-48 

44 

. 057 

— . 04i7 

1.307 

-.55 

54 

, 025 

-. 052 

1. 156 

-68 

39 

.008 

-.036 

1.098 

-80 

17 I 

. (KXl 

-. 029 

1,081 

-91 

59 

. 010 

-. 013 

1. 120 

-102 

5>3 

. 019 

. 000 

1.211 

-114 

55 

.036 

-. Oil 

1.370 

-128 

14 

.044 

-.040 

1.769 

-144 

16 

.020 

-.067 

2. 368 

-154 

50 

.000 

-.027 

19.83 

-179 

24 


TABLE II 


*(¥>>-0.1 sin (^-45°) ^o=0.10 /»=*(») -0.0657 = S'* 47' 


UPPER SURFACE 


ifi 1 

• 

0 


1 

dc 

1 

cosA 4' 

sinA ^ 

Degree.s 

Rartians 

Ra<lians 

Degrees 

(10 

0 

0.0000 

-0. 0707 

0. 0707 

4 3 

0. 1707 

0. 0707 

0. 0707 

1.0146 ^ 

0. 171.5 


.0873 

-.0643 

. 1516 

8 41 

. 17W) 

. 0766 

.0643 

1. 01.56 

. 17,5 

10 

. 1745 

-. 0574 

. 2.319 

13 17 

. 1819 

.0819 

.0574 

1.0166 

. 1826 


2618 

-, 0500 

.3118 

17 52 

. 186,6 

. 0866 

.0500 

1. 0175 

. 1877 

20 

. 3491 

-.042:1 ' 

.3914 

22 26 

. 1906 

.0906 

.0423 

1.0182 

. 1918 

25 

. 4.363 

-. 0.342 I 

. 4705 

26 57 

. 1940 

,0940 

.0:U2 

1.0189 

. 1952 

30 

. 5236 

0259 

. 549.5 

31 29 

. 1966 

,0966 

.02.59 

1.0194 

. 1979 


. 6109 

-.0174 

.628:1 

36 0 

, 1985 

.0985 

.0174 

1.0198 

. 1998 


78,54 

.0000 

. 7854 

45 0 

.2000 

. 1000 

.0000 

1.0201 

. 2013 

55 

. 9599 

.0174 

. 9425 

54 0 

. 1985 

.0985 

-.0174 

1.0198 

. 1998 

70 

1, 2217 

,(rt23 

1. 1794 

67 .35 

. 1906 

, 0906 

-.0423 

1.0182 

. 1918 

80 

1. 3963 

. 0574 

1. 3389 

76 43 

. 1819 

,0819 

-. 0574 

1.0166 

. 1826 

90 

1. 5708 

, 0707 

1. .5001 

85 57 

. 1707 

,0707 

-.0707 

1.0146 

. 1715 

100 

1. 74.53 

. 0HI9 

1. 6634 

95 18 

, 1574 

.'0574 

-.0819 

1.0124 

. 1581 

no 

1.9199 

, 0906 

1.829,3 

104 49 

. 1423 

. 042:i 

-.0906 

1.0101 

. 1428 

125 

2, 1817 

, 0985 

2. 08:12 

119 22 

. 1174 

i .0174 

-.0985 

1. 0069 

. 1177 

135 

2, 3.5<'>2 

, 1000 

2. 2.562 

129 16 

. 1000 

.0(X>0 

-. 1000 

1. 0050 

. 1002 

150 

2. 6180 

. 0966 

2.5214 

144 28 

.0741 

-.0259 

-.’0966 

1. 0028 

.0742 

160 

2. 7925 

.0906, 

2. 7019 

154 48 

. 0577 

-.0423 

-.0906 

1. 0017 

. 0577 

170 

2 9671 

. 0819 

j 2. 88,52 

165 19 

.0426 

-.0.574 

-.0819 

1.0009 

, 0426 

180 

3. 1416 

. 0707 

:i. 0709 

175 .57 

i 

. 0293 

-.0707 

-.0707 

1.0004 

.0293 


0. 9975 
.9885 
.9733 
. 9518 
.9243 
.8914 
.8,528 
.8090 
.7071 
. 5878 
. 3813 
. 2298 
.0706 
-.0924 
-. 2557 
-.4904 
-. 6329 
-. 8138 
-.9048 
-. 9673 
-. 9975 


0. 0703 
. 1510 
.2298 
.3068 
. 3816 
. 4532 
.5223 
. 5878 
. 7071 
.8090 
. 9244 
. 9733 
. 9975 
.9957 
.9,)68 
.8715 
. 7742 
, 5812 
. 4258 
. 2538 
. 0706 


1.0121 
1.0039 
. 9895 
. 9685 
.9411 
.9082 
.8693 
, 8250 
. 7213 
.5994 
. 3882 
.2336 
.0716 
-.0935 
-.2.583 
-. 4938 
-.6361 
-.8161 
-.9063 
-.9682 
-.9979 


0.0121 
.0268 
.0420 
. 0576 
. 0732 
. 0885 
. 1034 
.1174 
. 14-23 
. 1616 
. 1773 
. 1777 
. 1711 
. 1574 
. 1381 
. 1026 
.0776 
.0431 
.0246 
,0108 
.0021 


0 

0. 0(XX) 

-0. 0707 

0. 0707 

4 3 

-5 

-.0873 

-.0766 

-.0107 

-0 37 

-10 

174,5 

-.0819 

-. 0926 

-5 18 

-15 

-.2618 

-.08693 

17,52 

-10 2 

-20 

-. .3491 

-. 0906 

-. 2.585 

-14 49 

-25 

4.363 

-.0940 

342-1 

-19 :i7 

-30 

-. 5236 

-.09W3 

-.4270 

-24 28 

-35 

6109 

-. ai85 

5124 

-29 21 

-45 

78.54 

-. 1000 

-.68,54 

-39 16 

-.55 

-.9.599 

-.0985 

-.8614 

-49 21 

-70 

-1, 2217 

-.0906 

-1. 1311 

-64 48 

-80 

-1.3963 

-.0819 

-1.3144 

-7.5 19 

-90 

-1.5708 

-.0707 

-1..5001 

-85 57 

-100 

-1.74.53 

-.0574 

-1.6879 

-96 43 

-no 

-1.9199 

-. 0423 

-1.8776 

-107 3 

-125 

-2. 1817 

-.0174 

-2. 1643 

-124 1 

-135 

-2. 3.562 

.0000 

-2. 3,562 

-13.5 0 

-150 

-2.6180 

. 0259 

-2. 6439 

-151 29 

-160 

-2.7925 

. fM23 

-2.8348 

-162 25 

-170 

-2. 9671 

. 0574 

-3. 024.5 

-173 18 

-180 

-3. 1416 

.0707 

-3.2123 

-184 3 


LOWER SURFACE 


0. 1707 

0. 0707 

.0707 

1.0146 

0. 1715 

. 1643 

.0643 

. 0766 

1.0135 

. 16.50 

. 1574 

. 0574 

.0819 

1.0124 

. 1581 

. 1500 

,0,500 

.0866 

1.0113 

. 1506 

. 1423 


.0906 

1.0101 

.1428 

. 1342 

. 0: i 42 

.0940 

1.0090 

. 1346 

. 1259 

. 02.59 

.0966 

1,0079 

. 1262 

.1174 

.0174 

.0985 

1.0069 

. 1177 

. 1000 

,00( K ) 

. 1000 

1.0 a 50 

. 1002 

.0826 

-.0174 

.0985 

1.0034 

.0827 

. 0577 

-.0423 

.0906 

1.0017 

. 0.577 

,0426 

0.574 

! .0819 

1.0009 

.0426 

.0293 

-.0707 

.0707 
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6. 3941 
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4.0960 
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2. 8421 
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